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We study the local height probabilities in a composite height model, derived from the restricted 
solid-on-solid model introduced by Andrews, Baxter and Forrester, and their connection with confor- 
mal field theory characters. The obtained conformal field theories also describe the critical behavior 
of the model at two different critical points. In addition, at criticality, the model is equivalent to 
a one-dimensional chain of anyons, subject to competing two- and three-body interactions. The 
anyonic-chain interpretation provided the original motivation to introduce the composite height 
model, and by obtaining the critical behaviour of the composite height model, the critical behaviour 
of the anyonic chains is established as well. Depending on the overall sign of the hamiltonian, this 
critical behaviour is described by a diagonal coset-model, generalizing the minimal models for one 
sign, and by Fateev-Zamolodchikov parafermions for the other. 

I. INTRODUCTION 

Ever since the advance of conformal field theory in the seminal paper by Belavin, Polyakov and Zamolodchikov [IJ , 
it has played an extremely important role in the study of critical behaviour in two-dimensional statistical mechanics 
models, and one-dimensional quantum systems alike. Not only were the foundations of conformal field theory (CFT) 
laid out in £Q, in addition, an infinite series of conformal field theories were introduced, the so-called minimal models. 
These CFTs describe an infinite number of possible critical points. Roughly at the same time, Andrews, Baxter and 
Forrester [2J studied a generalization of the eight-vertex model, in which the degrees of freedom are heights, living 
on the square lattice. These heights can take a finite set of r — 1 values, where r is a parameter characterizing the 
model. Because of this constraint, these models are also called 'restricted-solid-on-solid' (RSOS), or simply 'height' 
models. These models were shown to exhibit various gapped phases, separated by critical points. Shortly afterwards 
Huse [3] realized that the critical points found by Andrews, Baxter and Forrester are described by the family of 
unitary minimal models obtained by Friedan, Qiu and Shenker around the same time [3]. One can certainly say that 
conformal field theory took a flying start! Here we study aspects of a similar connection between critical points in 
composite height models, as well as in anyonic quantum chains, and CFTs. 

Specifically, we study the local height probabilities (LHPs), to be defined below, of generalized RSOS models [5] 
inspired by an anyonic quantum chain with competing two- and three-body interactions, which was introduced in [BJ. 
Originally [7], the anyonic chains were motivated as simple models for interacting anyons in topological phases, hosting 
the anyons constituting the chain as their elementary excitations. Indeed, the local degrees of freedom of the anyon 
chain related to the height model we study here are non-abelian su(2)k anyons and — exactly as for a Heisenberg 
chain of su(2) spins — the anyon hamiltonian assigns an energy cost depending on to which representation the two 
or three neighboring anyons are 'fused'. In contrast to spin chains, however, only the k + 1 integrable representations 
of su(2)k can appear and, crucially, the Hilbert space of the chain is not a local tensor product of the single anyon 
degrees of freedom. Previous studies have revealed that the anyonic chains have very rich phase diagrams, even 
richer than in the original spin models, with novel phases protected or broken by 'topological symmetries' [BJ [7]. 
Here we study two integrable critical points of the anyon model with competing two- and three-body interactions, 
identified in [5]. These integrable points of the anyon chain are directly obtained from an 'anyonic' representation 
of the Tcmperley-Lieb algebra and are equivalent with the critical points of a classical, integrable 'composite' RSOS 
height model also introduced in [5]. This mapping is crucial, since the non-local form of the Hibert space of the anyon 
chain reduces the utilizability of Bethe ansatz techniques (but see [H [9]). 

The RSOS height model is defined on a square lattice with the heights k being local degrees of freedom on each 
vertex i, subjected to the constraints 1 < li < r — 1 with \li — lj\ = 1 for adjacent vertices i, j, and r — k + 2 for su(2)k 
anyons, and is a composite of the original integrable RSOS models of Andrews, Baxter and Forrester (ABF) [2 . The 
various critical points of the ABF model, separating the ordered phases, were subsequently identified in [3], and shown 
to provide realizations of the minimal models studied in the seminal papers [TJ 0] . This type of multi-critical behavior 
also applies to the composite model. The study of the critical, as well as off-critical behavior of the model goes via 
the so-called local height probabilities, which are the probabilities for a central site to have a certain height, given the 
boundary conditions. As with the original ABF model, the off-critical LHPs of the height model exhibit properties 
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related to conformal field theory characters and can be calculated exactly with the corner transfer matrix method 
[T01 - H2] . This off-critical CFT structure of the LHPs is governed by the same theory that describes the critical points 
of the lattice model; the former arises from integrable perturbations of the latter. In particular, for a finite lattice 
size, the LHPs are composed of finitized forms of CFT characters from which one can obtain the characters by taking 
the thermodynamic limit. 

Thus the two integrable points of the anyon chain are related to two different regimes of the height model, regimes 
II and III in the notation of [H[S], not just the critical points to which they terminate. Therefore LHPs allow one to 
determine the (extended) critical behavior of the anyon chain once the off-critical CFT has been identified, which is 
the aim of the current paper. 

The central objects of interest in our paper are fermionic generating functions, quantities usually called universal 
chiral partition functions (UCPFs) [13] , or fundamental fermionic forms p31 [T5] , which we want to reproduce the 
LHPs in closed form and relate to finitized forms of CFT characters. 

As with the first proofs of the connection with the ABF model and the M.{r — 1, r) minimal models, our strategy of 
proof is based on the recurrence properties of the polynomial UCPFs and LHPs, this strategy is sometimes referred to 
in the literature as Schur's method. In the thermodynamic limit, our formulas give fermionic characters of the coset 
conformal field theory s "( 2 h x ^( 2 )i x s"( 2 )r-4 m ^ ne regime III and Z r -2 parafermions in the regime II, corresponding 
to the two integrable points under study. Based on numerical checks [S[ [S] , these CFTs were earlier identified as the 
critical behavior of the anyonic chain. The parafermion theory and characters are well known and appear also in the 
ABF models; we recover the fermionic forms of these characters [T5H20|. The fermionic characters we obtain for the 
coset theory are new to the best of our knowledge. 

Finally, the fermionic forms of the LHPs of the composite model studied here open an arena of q-identities related 
to the coset theories sn ( 2 )i x ^( 2 )i x s"( 2 )r-4 _ T ms } s exactly like the g-identity and CFT character results obtained for 

the minimal models from the ABF-type models [21 HH [HI [17l [HI I2TH27] and provides another motivation to study 
the composite height model. The UCPFs that we obtain in the regime III for the coset theory are characterized by 
the fact that they have two 'real' fermions and r — 5 'pseudo' particles, instead of just one 'real' fermion appearing in 
the ABF type expressions [221[23H2S]- Given that one would be able to obtain the bosonic forms of the characters, 
giving interesting bose-fermi type identies, one could possibly also obtain new types of Rogers-Ramanujan and q- 
identities 2, 9] Q31 [T51 US US] . As already mentioned, the dual finitized characters obtained for the regime II are Z r _ 2 
parafermions, as in the ABF case, with all fermions 'real'. In fact, this type of behavior seems to be rather generic, 



as our calculations in sec. VI B suggest. 

This paper is organized as follows. In scction[H]we briefly introduce and recollect the composite height model from 
Ref. [5] and set out the stage and notation for the various quantities related to the LHPs. As in the height model of 
ABF, the composite model is parametrized by an integer r that sets the maximum of the height variables, with r > 5 
for the anyonic chains. In the paper [5], a fermionic form for the central quantity X m in the LHPs for the simplest 
case r = 5 was introduced, where m is the lattice size of the height model, based on numerical checks. In section [Hi] 
we prove this equality analytically and then proceed for the r = 6 case in sec. |IV[ In both cases, the proof is obtained 
using the recurrence properties of the LHPs and the fermionic generating functions related to those with respect to 
the lattice size, exactly as in the original ABF model. In section [VJ based on the structure for r — 5,6, we give the 
general form of our fermionic UPCFs related to the LHPs, a conjecture we claim valid for any r and supported by 
numerical checks for r > 6 and correct central charges. While we could apply the same the strategy of proof for bigger 
r, the computations become quickly cumbersome and not very illuminating. In sections VI and |VIl[ we finally study 



the thermodynamic limit of the UPCFs and give the explicit connection with CFT characters, respectively. We end 
by discussing our results and giving some future directions of study. 

II. LOCAL HEIGHT PROBABILITIES IN A COMPOSITE HEIGHT MODEL 

The height model we study in this paper is most easily explained in terms of the original height model introduced 
by Andrews, Baxter and Forrester (ABF) [2J. We will first introduce this model, and subsequently explain how the 
composite height model can be constructed from it. 

A. Definition of the height models 

The ABF model consists of heights which are assigned to the vertices of the square lattice. The heights I can take 
the values I = 1,2, ... ,r — 1, where r is an integer satisfying r > 3. Different values of r correspond to different 
models; the case r = 5 is equivalent to the hard hexagon model, as explained in [2]. Heights at neighbouring vertices 
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have to satisfy the constraint that they differ by one. This constraint leads to six different types of plaquettes on the 
square lattice, as depicted in Fig.JT] The two plaquettes labeled by ai (and similarly for /3j) will be assigned the same 









FIG. 1: The six different type of plaquettes occurring in the ABF model. 



weight, so that one obtains an isotropic model. ABF showed that this model can be solved for a two-parameter family 
of weights. We will denote these parameters by p and u. The first parameter — 1 < p < 1 resembles a temperature, 
and drives a phase transition at p = 0. The parameter u is related to the 'anisotropy' of the lattice, and is the variable 
appearing in the Yang-Baxter equation below. The behaviour of the model will not depend on the magnitude of u, 
only its sign. For a description of the various phases of the ABF model, we refer to the original paper [2] , and the 
paper by Huse [3J, who studied the connection between the critical points of the model and conformal field theory. 
We will come back to the various phases of the composite height model after we explained how the model can be 
solved. 

The weights for which the ABF model can be solved explicitly are given in terms of elliptic functions. To specify 
them, we introduce the following notation. First, p will be related to the modulus m 2 of the theta functions via 

_ K'(tti) 

p = e n A " (m > , where K(m) is the complete elliptic integral of the first kind and K'(m) = K(l — m). The parameter 
r enters the weights via r] — K ^ , while the values of the heights I enter as W[ — 2r\l. 

Introducing the elliptic functions H(u) — 6\ ( 2 K(m) > P) an< ^ = ®^ 2K{m) we define h(u) — H(u)Q(u), where 
we have suppressed the dependence on p (or m), as is customary. The 9i are the Jacobi theta functions. Explicitly, 
one finds the following expression for h(u), 



h(u) = 2pi sin(^) 2 J3 "cos(^) + p 2 ")(l -p 2 ") 2 
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We can now introduce the two-parameter family of plaquette weights as follows 

h(u) (h(wi~ 1 )h(wi +1 )) [ 
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As a first step in solving their model, ABF noted that the weights ^ satisfy the Yang-Baxter equation, which 
implies that the row-to-row transfer matrices for different values of u commute with each other. We will explain in a 
bit more detail how the model was solved below, after we introduced the composite height model considered in [5]. 

Inspired by the the work of Ikhlef et al. [551 on loop models, a composite model was constructed in the following 
way. First, we denote the weights associated with a plaquette by W(u; l\,h,h^ U) as depicted in Fig. [2j 

By taking four of these (2 x 2) plaquettes, one can from a composite (3 x 3) plaquette, with the weight (see Fig. [3]) 

W(u; h,h,hM, h, h, h, Is) = W ( u ' h, h, I, k)W(u + K- l 2 ,l 3 , l 4 , l)W(u; I, l 4 , l 5 ,h)W{u - K; l 8 ,l, l 6 ,l 7 ) ■ (3) 



We note that the parameter u of two of the sub-plaquettes has been shifted. Without this shift, the composite model 
would be equivalent to the original model. It is a straightforward exercise to show that the plaquette weights W of 
the composite model satisfy the Yang-Baxter equation, by only making use of the fact that the plaquette weights W 
of the original model satisfy the Yang-Baxter equation, namely 



W(u; l u l, l 5 , l 6 )W{u + v; l 2 , l 3 , 1, h)W(v; l 3 , i 4 , Z 5 , 1) = ^ W(v; l 2 , 1, k, h)W(u+v; I, l i} h,l 6 )W(u; l 2 ,l 3 , l 4 , 1) . (4) 
i i 
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FIG. 2: The plaquettes of the ABF model. 



B. A glimpse on the Corner Transfer Matrix method 

We now briefly discuss how these height models can be solved. One makes use of the Corner Transfer Matrix 
(CTM) method, which was explained in detail in chapter 13 of Baxter's book [3D] (see also [3T] for a recent account). 
The application of the CTM method to solve the ABF model is detailed in [2], while [5] deals with the composite 
model. 

The key objects in the CTM method are four corner transfer matrices. In contrast to the row-to-row transfer 
matrix, the corner transfer matrices do not add merely one row to the lattice, but instead an entire quadrant, or 
corner. If we denote the CTMs of the four different corners by A, B, C and D, the partition function of the model is 
given by Z — Tr(ABCD). Using the CTM method, one can calculate a quantity called the local height probability 
(LHP). In particular, P a denotes the probability for the central height, say l\ (at the corner of the four CTMs) to 
take the value a. This probability can be written as 



P„ 



Tr(S a ABCD) 



(5) 



where S a is the diagonal matrix, with diagonal entries 1 if the central height li = a, and zero otherwise. 

To make this discussion a bit more explicit, we display the CTM A explicitly in Fig. [4] The rows and columns of 
A are labeled by (I1J2, ■ ■ ■ lm) and (l[ = h,l' 2 , ■ ■ ■ , l' m ) respectively. The central height li is fixed to be a, while the 
heights l m +i, l' m+ i, l m +2, l' m +2i etc., at the boundary are fixed to ground state values of the model. By analyzing the 
weights of the model, one can show [5] that the ground states are in fact diagonal, in this case along the SW-NE 
direction, and are fixed by the boundary conditions (6, c, d, e) as indicated in the figure. The different ground state 
'patterns' (&, c, d, e) are discussed in detail in section |TlD[ following [5]. The CTM method allows one to calculate the 
local height probabilities P ai and these depend on the boundary conditions (b,c,d,e). The matrices B, C and D are 
obtained in a similar way as A, by subsequent rotations of the diagram over 7r/2. From the definition of A, it is clear 
that the 'size' of the quadrant m, equal to twice the number of the top-row plaquettes plus one, has to be odd. 

We do not explain the calculation of the LHPs in full detail, but merely state the main ingredients (following 0E]) 
of this calculation and of course the result, which is an expression for the height probabilities. These local height 
probabilities are the starting point for the current paper, and the goal is to prove that the LHPs are equal to (Unitized) 
characters in conformal field theory (CFT). 

The first essential ingredient of the CTM method is to make use of the Yang-Baxter equation, to show that the 
CTMs A etc. can be written in a special, diagonal form, see [30] for the details. One starts by formally equating the 
product lim m ^oo B(u)C(v), which covers half of the lattice, to the limit lim„_ i . 00 T(u, v) n using the inhomogenous 
row-to-row transfer matrix T(u,v), the latter limit also then covering half of the lattice with the anisotropics u,v 




E 




FIG. 3: The plaquette weights of the composite model, where the symbols ± denote the corresponding shifts of it in j3|. 
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FIG. 4: The CTM A of the composite model, with the boundary conditions (b, c, d, e) as specified in the text. 



in the two quadrants [30] . The Yang-Baxter equation can be shown to ensure that B{u)C(v) only depends on the 
difference u — v. In the end, one obtains the following form for the CTMs 

A(u) = Q l M l e- uU Q^ 1 (6) 
B(u) = Q 2 M 2 e uU Q^ 
C(u) = Q 3 M 3 e' un Q^ 
D(u) = Q 4 M 4 e un Q^ , 

where H, Qi and Mi (with i = 1,2, 3, 4) do not depend on u, commute with the matrices S a , and in addition H and 
Mi are diagonal. 

To relate these diagonal forms of the corner transfer matrices to the height probabilities, one has to calculate the 
form of the CTMs for particular values of u. In the case of the composite model, one has to use an identity relating 
particular sums of products of elliptic functions to a single product. The details can be found in [5] but we summarize 
the results here. By the periodicity properties of the elliptic weights, one only needs to consider u in the two following 
regimes 

X>i : < u < 2-q + K = (2 + r)r), V 2 : 2?? - K = (2 - r)-q < u < 0. (7) 
One can show that, up to scalar multiples, 

A(0) = Q^hQz 1 = 1, (8) 

which allows one to write 

A{u) = Q 2 e- uH Qz\ (9) 
so the diagonal form of A(u) is equal to an exponential. The height probability P a is then given by 

P n (h,c,d,e) = ^ MlM2M3M f ) . (10) 
y ' ' ' ' Tr(MiM 2 M 3 M 4 ) V ' 

and will depend in addition on the boundary conditions (b, c, d, e). The product M\M 2 M 3 M 4 can be computed by 
considering different limits of the corner transfer matrices. First, in the limit u — > in the domain T>\ and up to 
irrelevant scalar factors, one has 

A(0) = C(0) = 1, (11) 

and secondly, when u —¥ (2 + r)rj, one has 

B(u = (2 + r)rj) = D(u = (2 + r)r?) = V x , (12) 

where 



(Vl\v - y/h{2t]h)S(l,l'). 



(13) 
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Therefore, 

A(0)B{{2 + r)r))C(0)D((2 + r)rj) = M 1 M 2 M 3 ,M i e^ 2+r ^ = V?. (14) 
Similarly in the domain T>2, where the weights effectively only change their signs, 

A(0)B({2 - r)rj)C(0)D((2 - r)r)) = MyM^M^M^-^ = V? . (15) 
These give the height probability as 

{ ' Tr^e- 2 *^) \2-r, u e V 2 V ' 

The final step in the calculation is determining the diagonal form of the CTMs. To do this, one first employs the 
'conjugate modulus transformation', which gives an expansion of the weights around p = 1, instead of p — in the 
original formulation of the weights. The result of this calculation is that the CTMs are diagonal in the limit p —> 1. 
In calculating the matrix elements, the first observation is that the elliptic weights of the model are quasi-periodic in 
u with the period 2iK' , so the elements of % are integer multiplets of tt/K' , 

«,,„ = !5H5, (17 , 

where N(l) is an integer function. In particular, A takes the form Ay = {e~ u ^ i )i t v — gj^ w^'Si^i , where 1 = 

(li,...,l m ) and 1' = . . . ,l' m ) label the rows and columns of A; w — e~ 27r x T , and gi x = w ( ibf - . Finally 
= N(l)/2 is given by 

m + l 

M = E i ( lhj+3 4 hi ~ A + ^-a^t 2j+1 M )+3 Si 2j ,i 2j+ ^ . (is) 
j=l \ / 

Having found the diagonal form of A in the limit p — > 1, one uses the last essential ingredient of the method, to find 
the diagonal form for all (positive) p. The function <fi(l) takes integer or half-integer values. Because the weights of 
the model depend continuously on p, it is reasonable to assume that W does not change discontinuously with p. This 
in turn implies that the function 0(1) is in fact independent of p, and the diagonal form of A which was determined 
for p = 1, is in fact valid for < p < 1. With this diagonal form for A, one can give an explicit expression for the 
local height probabilities P a - 

C. The local height probabilities and the function X m (a;b,c,d,e;q) 

The local height probabilities can finally be written in the following from 

P a (b,c,d,e) = S^VaXmi^b, c,d, e;x*). (19) 
Va = x^-^^/^Eix^x^ (20) 
S = ^v a X m (a;b, c, d, e;x l ) (21) 



= p -in V /K' = -i^K/K' 



(22) 



with boundary conditions l\ — a and l m +i, Un+2, • ■ ■ = b,c,d,e and m is the lattice size. The variables p and x both 
lie in the range < x,p < 1, but when p — > 0, we have x — > 1, and vice versa. The function E(z,x) is the triple 
product 

oo 

E(z,x)= Y[(l - x"- 1 z)(l - x n z- 1 )(l - x n ). (23) 

n=l 

The height probabilities take different forms dependent on the parameter u, which enter the expressions for the 



probabilities via t, as indicated in (16). In the case that u > 0, which we call 'regime III', following the notation in 
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[2J, t = r + 2, i.e. t is greater than zero. In the case u < 0, 'regime IF, we have that t — 2 — r, i.e. t is less than zero. 
For more details on the regimes, we refer to [2j[l]. 
The function X m (a; b, c, d, e; q) is defined as 

X m (a;b,c 1 d 1 e;q)= ]T (24) 

1— (a, 02 )■ ■ - 7 lm )fc)C,d,e) 



where 0(1) is the function defined in ( 18 1. The boundary conditions were l\ — a and l m +i — b, . . . , l m +A — e, and the 
lattice size m is odd in the composite model. The heights are diagonal in the limit p — > 1, and from the definition 
one can see that 

X m (a; b, c, d, e; q) = g =^( ] * = i= 5L +*-i..*«..*M) X m _ 2 (a; 6 - 2, 6 - 1,6, c; g) 

+ g 2 * i (M= £l +^-i,o5c, e 5b,< J )x m _ 2 (a; 6, 6 - 1, b, c; g) (25) 

+q ^{^^+8 b+1 _ c S^S b<d ) Xm _^ a . b ^ b+h 6j c . g ) 

+g^(^ ± ^ +5b + 1 ^- A ")x m _ 2 (a; 6 + 2, 6 + 1, 6, c; g) 



and the states appearing on the right hand side (RHS) only depend on b, c. Again, by assuming continuity in p, this 
recursion relation is valid for all p. Also, due to the symmetries of the plaquette weights of the model, X m satisfies 
015] 

X m (r - a;r - b,r - c,r - e;q) = X m (a; b, c, d, e; g). (26) 



D. Phases of the composite height model 

From the expressions of the local height probabilities (or better, the partition function), we can extract the phase 
diagram of the model. Here, we will concentrate on the case p > 0, for which the expressions for the FHPs of the 
previous section are valid. 

In this section, we give the ground states in the gapped region < p < 1. These ground states will play an 
important role in making the connection between the FHP for p — and the conformal field theory characters. 

We start by considering the case u > 0, i.e. regime III. The ground states are those configurations which contribute 
maximally to the partition function. As was the case for the FHP, the dependence on the regime is via the parameter 
£, which is positive in regime III. This in turn implies that to find the ground state configurations, the function 0(1) 
has to be minimized (see J3 for more details). The first term in 0(1) vanishes when hj+3 = hj-i, which is a necessary 
condition in a ground state. We recall that neighbouring heights have to differ by one. The first way in which second 
term in 0(1) also vanishes is when hj = hj+2 and hj+i = hj-i ± 2. These type of ground state patterns will be 
denoted by Gf (when hj+i = hj-i + 2) and Gj (when hj+i = hj-i — 2). The second term in 0(1) also vanishes for 
hj+i — hj-i an d hj — hj—i + 1 = hj+2 + 2 (these ground states are denoted by Gj), or for hj — hj-i — 1 = hj+2 — 2 
(these ground states are denoted by G^)- This exhausts the possible ground state patterns for u > 0. These ground 
state patterns are depicted, together with the other possible patterns to be discussed below, in figure [5j 

Turning our attention to the case u < or regime II, we have that the parameter t is negative, which implies 
that the ground states maximize the function 0(1) (see [5]). We start by noting that the first term in the sum in 
0(1) can not be one for all values of j, because the heights take their values in the finite range I = 1,2, ...r — 1. 
The second term in the sum can however always be one, which is thus the case for the ground states. We find that 
hj—i = hj+i = hj+3 and hj — hj+2 = hj-i ± 1 are the necessary conditions. These u < ground state patterns are 
denoted by by G 3 . 

Before discussing the transition point p — 0, we first mention that, if one considers four consecutive heights, there 
are only two patterns left, which are not part of a ground state pattern. These are (6, b + 1, b + 1, b + 3), which we 
will denote by NGS + and (b, b — 1, b — 2, b — 3), denoted by NGS~ . These patterns will play a role in the study of the 
local height probabilities in connection with the critical point at p = 0. 

At the point p — 0, we find that all configurations contribute to the partition function, which means that one has 
to consider the local height probabilities in their entirety. In [5J , it was observed that the local height probabilities are 
related to (finitized) characters of certain conformal field theories. This behaviour, as alluded to in the introduction, 
has been observed in various other cases as well [T0U12) . not in the least for the ABF model. 

Fet us be a bit more precise about the connection between the expressions for the FHPs and conformal field theory. 
It turns out that the functions X m (a; b, c, d, e; q) appearing in the expression for the FHP P a (b, c, d, e) correspond 
to a character of a conformal field theory, if the boundary condition (b, c, d, e) is part of a ground state pattern we 



^,6 + 3 NGS+ = (6,6 + 1,6 + 2,6 + 3) *6+l = (6,6- 1,6,6+ 1) 

6 + 2 x ? . h ' \h<r-\ 

G+ = (6,6+1,6 + 2,6 + 1) \ ^ \ G 3 - = (6,6-1,6,6-1) 



6+1 6+1 i& = 



6-1 6-1 Jt = r-1 



& ^ ^ b S G+ = (6, 6 +1,6, 6+1) \ ^ Gr = (6,6- 1,6-2,6-1) 

^6-1 = (6,6+1,6,6- 1) ^6-3 NGS" = (6, 6 - 1, 6 - 2, 6 - 3) 

FIG. 5: The various height configurations that appear in the model. The patterns G^ and G^ are ground states for u > 
and the patterns G^ for u < 0. The patterns NGS* are not ground state patterns but appear in the recursion of X m . The 
meaning of the arrows is explained in the main text. 



discussed above. The conformal field theory is the one describing the critical behaviour at the phase transition to the 
phase exhibiting the ground state pattern under consideration. 

In particular, in [5], an explicit expression for the function X m (a; b, c, d, e; q) for r — 5 was conjectured, which 
equals the Unitized characters of a particular CFT. In the following, we prove this conjecture, thereby establishing 
the connection between the model for r = 5, and the CFT. In [5], the connection with the Gepner parafcrmions 
associated with su(3)2 was made. This parafermionic coset sw(3)2/(u(l)4 x 14(1)12) is equivalent with the diagonal 
coset su(2)i x su(2)i x su{2)\/ su(2)^, which is only one member of an infinite series of equivalences, which starts with 
the equivalence of the Z 2 parafermions, 57/(2)2/^(1)4 and the first minimal model, i.e. the Ising model. In addition, 
we provide an explicit form for the functions X m (a;b,c 1 d 1 e;q) for arbitrary r, and prove the result also for r = 6. 
We argue that these functions are the Unitized characters of a set of coset models, similar to the 'minimal models' 
describing the u > critical point of the ABF model. 

Before we start the discussion of the general properties of the functions X m (a; 6, c, d, e; q) in the next subsection, we 
note that the number of independent height probabilities, given the boundary conditions which correspond to ground 
state patterns, is (r — l)(r — 3) in the regime u > 0, and (r — l)(r — 2)/2 in the regime u < 0. Here, the reflection 
symmetry (26 1 was already taken into account to reduce the number of independent functions X m (a; 6, c, d, e; q). 



E. Recursion relations for X m (a;b, c, d, e; q) 

We continue by describing some general properties of the functions X m (a] b, c, d, e; q), before we deal more explicitly 
with the cases r — 5 and r = 6 in the following sections, where we give an explicit expression for these functions, and 
prove that they are equivalent to the functions X m , by showing that they obey the same recursion relations, and have 
identical boundary conditions. 

The possible boundary conditions (b, c, d, e) for the LHP are only constrained by the fact that neighbouring heights 
have to differ by one, and are all given in figure [5] 

The different states labelled Gf,Gf,Gf and NGS ± appear in the recursion for X m as described below. 

Although in establishing the connection between the LHPs and conformal field theory, we are mainly interested in 
the boundary conditions corresponding to ground states, the recursion relations force us to consider the non ground 
state patterns as well, because they are generated by the recursion relations automatically. 

Due to the relation ( 26 ) , we need to consider only half of the patterns in figure [5j We will focus on the patterns 



with an increasing second height. Starting from an G± boundary condition, relevant for u > 0, one finds 

X m (a; b,b + l,b + 2,b + l;q) = q^ X m _ 2 (a; b - 2, b - 1, 6, b + 1; q) +q^ X m _ 2 (a; b, b - 1, b, b + 1; q) (27) 



NGS+ Gz 

+ X m _ 2 (a; 6, b + 1, 6, b + 1; q) +X TO _ 2 (a; b + 2, b + 1, b, b + 1; q) . 

v ' v ' 

G t G l 

Here the first term vanishes for b < 2 and the second in the case 6=1. The third term is a ground state pattern for 
u < and the first term is a non-ground state pattern. 
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For the ground states in G^ , we get 

X m {a- 6,6+1,6,6— = X m _ 2 (a; b - 2, b - 1, 6, 6 + 1; q) + X TO _ 2 (a; 6, 6 - 1, 6, 6 + 1; g) (28) 

V v ' v v ' 

NGS+ Q- 

+q^ X m _ 2 (a; 6, 6 + 1, 6, 6 + 1; q) +q^ X m _ 2 (a; 6 + 2, 6 + 1, 6, 6 + 1; g) . 



Gj 

Again, some terms do not necessarily contribute and the first term is not a ground state pattern. 
The set of patterns relevant for the u < LHP lead to a recursion of the form 

X m (a; 6, 6 + 1, 6, 6 + 1; q) = q^r 1 AT m _ 2 (a; 6 - 2, 6 - 1, 6, 6 + 1; q) +q m f 1 X m - 2 (a; 6, 6 - 1, 6, 6 + 1; q) (29) 



NGS+ G 2 

+q m t 1 X m _ 2 (a, 6, 6 + 1, 6, 6 + 1; q) +q n ^ X m _ 2 (a; 6 + 2, 6 + 1, 6, 6 + 1; g) 



gJ g ± 
Finally, the non-ground state pattern NGS + satisfies the the recursion 

X m (a; 6, 6 + 1, 6 + 2, 6 + 3; q) = X m _ 2 (a; 6 - 2, 6 - 1, 6, 6 + 1; q) +q m ^ 1 X m _ 2 (a; 6, b - 1, 6, 6 + 1; g) (30) 



NGS+ 

X m _ 2 (a; 6, 6 + 1, 6, b + 1; g) X m _ 2 (a; 6 + 2, 6 + 1, 6, 6 + 1; g) . 

^ V V 

g 3 + cr 

It is clear that not only does the recursion relation generate non ground state patterns, it also mixes the different 
ground state patterns relevant for the two different regimes u > and u < 0. This implies that we will have to 
establish the recursion relations for all types of boundary conditions. 

To make our work easier, we start by establishing some relations between the different functions X m (a\ 6, c, d, e; g), 
which reduces the number of cases we have to check explicitly. These relations originate in the fact that the states 
appearing in the recursion only depend on 6, c. Using the definition of X m (a; b, c, d, e; q), one can see that changing 
the boundary height e — > e ± 2 as follows 

(6,6+1,6 + 2,6 + 1) ~* (6,6+1,6 + 2,6 + 3) 
(6,6 + 1,6,6 + 1) ~+ (6,6+1,6,6-1) 

leads to the relations 

X m (a;6,6+l,6 + 2,6 + 3;g) = q^ X m (a; 6, 6 + 1, 6 + 2, 6 + 1; g) (31) 
X m (a; 6, 6 + 1, 6, 6 + 1; g) = q^X m (a; 6, 6 + 1, 6, 6 - 1; g) , (32) 

where the second equation is only valid for 6 > 1. The first equation relates the the non ground state pattern NGS + 
to the u > ground state pattern Gf. The second equation relates the u < ground state pattern G3 to the u > 
ground state pattern G^~, for 6 > 1. The case 6=1 can be dealt with by relating the pattern to the u > pattern Gf 
by changing the boundary height d instead, namely (6, 6 + 1, 6 + 2, 6 + 1) (6, 6 + 1, 6, 6 + 1), as follows, (see fig. 

X m (a; 1, 2, 1, 2; q) = q^ X m (a; 1, 2, 3, 2; q). (33) 

We note that this relation only holds for 6=1, and can not be used to reduce the number of independent functions 
X m (a; 6, c, d, e; q) even further. In conclusion, we find that all the functions X m (a; 6, c, d, e; q) corresponding to non 
ground state pattens and u < ground state patterns can be related to u > ground state patters, and we are thus 
left with (r — l)(r — 3) independent functions, corresponding to, say, the G^~ and G 2 patterns. 



III. EXPLICIT EXPRESSIONS FOR r = 5 



We start our search for explicit expressions for the functions X m (a;b,c,d,e;q) with the case r — 5. For this case, 
an explicit functional expression was obtained in [5] based on numerical evidence, but the equivalence was not proven. 
We provide the proof in this section. 
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A. The function y(k;l2,h,l4;q) for r — 5 



For r = 5, the functions X m (a; b, c, (i, e) for the different boundary conditions are related to the functions [5] 



y{k;l 2 ,h,U;q) = 2^ q 

mi,T7l2>0 

1,...,4, Z 2 , ^3 = 1,2 and 



where L 



j(mj+m a - niimj- mi <5i 4) 3— m2#! 4 , 2) 



fc+rre 2 + ^i ai i+jj4 i 3 



2 

mi 



'fc+mi+<5j 2 ,i+<5; 3 , 2+^ 4 ,2 
2 " 
7772 



(34) 



777 




777 


77 




77 



is the g-binomial coefficient First, we define (q) m = YYjLiO-~q J ) 



for integer 777 > and (q)o = 1. We then define the (/-binomials, non-zero for integer m,n, as 

(g)m 



777 


."I 


77 













if < 77 < 777 integers, 
otherwise 



(35) 



The prime on the summation in (34) indicates constraints on the parities of 777 1 , 777 2 but with the above definition of 
the g-binomials, they are in fact superfluous and will be therefore left implicit. This, however, will not be the case for 
the constraints for r > 5. 

The relations between the X m and y{k) are schematically [Sj 



777 — 1 

X m (a;b,c, d,e;q) ~ y{ — - — ,l 2 ,l 3 ,l 4 ;q), 



(36) 



where the odd integer 777 is the lattice size of the composite height model and Z2, I3, I4 are determined by the configu- 
rations (a; b, c, d, e). 

Given this correspondence, the recursion for X m (a; b, c, d, e; q) implies a recursion schematically of the form 



y( m - ;l2,h,k;q) ~ <i 



or more conveniently in terms of the integer k = 



(™+l)/4 y( HLl _ 1; ^ l' 4;q) + ... + y(^-l 1; /»; ? ) 



m— 1 
2 ' 



la, h, U; «) ~ g -8 w(* - 1; «a, l' 3 , ii; «) + ••• + w(fc - 1; ia. 's. «) + 



(37) 



(38) 



for the functions y(k; la, I3, 14; g). From now on, we will display the dependence on the size of the lattice through the 
(integer valued) variable k instead of the odd integers 777. 

As in the original paper of ABF, we now set out to prove that the functions A 2 / £ +i(a; b, c, d, e; q) and y(k; l 2 , 13, 14; q) 
satisfy the same recursion relations and have identical boundary conditions and thus have to agree identically. This 
verifies the critical properties of the anyon model, since the functions y(k; la, I3, 14; q) are — conjecturally for general 



Unitized forms of CFT characters of the coset 



su(2)i X s u(2)i X su(2),._4 
su(2) r 



as in the original case of RSOS height 



probabilities and minimal models M(r — l,r) studied by ABF and by many authors in subsequent papers cited in 
the introduction. However, due to the composite nature of the height model, the corresponding recursions are more 
complicatedly related, as we will see. Unfortunately, we have not been able to obtain a functional form that would 
directly satisfy the recursion of X 2 k+i and need to proceed in a more oblique way in terms of the more general 
functions j/(fci, k 2 ; la, I3, 14; q), presented in the next subsection. The reason behind this is that the recursion for 
X 2 fe + i(a; b, c, d, e; q) gives a sum in terms of X 2 k-i, but all with different boundary conditions, while the recursion 
for y(k±, k 2 ; l 2 , 13, 14; q) leads to a sum of functions with the same values of the la, I3, I4, but with different values for 
h, k 2 . 

We first deal with the simplest case r = 5, corresponding to the diagonal coset su(2)i x su(2)i x su(2)i/su(2)3 (or, 
equivalently, the su(3)a/(u(l)4 x u(l)i a) G epner parafermions) in the regime u > and Z3 parafermions for u < 0, 
as initiated in the paper [5]. In section IV we prove the correspondence for r = 6 and give the general conjecture for 
arbitrary r in section [V] 



Recursion for y(k; I2, 13, q) 



The function y{k;l 2 ,lz,li;q) satisfies a recursion relation in k, based on the recursion for g-binomial coefficients 
which follows directly from the definition [32] 



m 




777 —1 




777 —1 




777 —1 




777 —1 


77 


= q n 


77 


+ 


77—1 




77 


+ q m ~ n 


77—1 



for 777 > 77 > 1 



(39) 
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In order to have a closed recursion for y(k; l 2 , h, h', q), we define a closely related function, which we will also denote 
by y and hope that no confusion arises, 



m 2 ,— mi m 2 — mi5i l<3 — m 2 5i 4 ,2) 



mi ,7712 >0 



"fei+m 2 +5i 3 ,i+<5( i 
2 

mi 



fc 2 +roi+5i„,i+(5j ,.2+5i, 2 



2 

m 2 



(40) 



Clearly y(k; l 2 , h, k] q) = y(k,k;l 2 ,l 3 ,l4;q). 

Then, using the latter recursion in (39) for kx or k 2 , leads to 



y(kx,k 2 ;h,h,h;q) = q 2 y(h - 2,fc 2 + i;l 2 ,h,h;q) + y(h - 2, k 2 ;l 2 ,l 3 ,U;q). 

k 2 + s l 3 ,2+ s l 2 ,l~ 1 

= q 5 y(h +l,k 2 - 2;l 2 ,l 3 ,k;q) + y(h,k 2 - 2;l 2 ,l 3 ,h;q). 



(41) 
(42) 



This is very similar to the recursion for a closely related function function Y(kx, fc 2 ; q) studied in Ref. (33j . Similarly, 
we can use the recursion in both k\ and k 2 to arrive at 



y(h, k 2 ;l 2 ,l 3 , U; q) = q s y{h -l,k 2 - 1, l 2 ,l 3 , U; q) + y(h - 2, k 2 - 2, l 2 , l 3 ,h; q) 

t l+ i l;,l+ i l 3 ,2~ 1 fc 2+ a i 3 ,l- 1 

+q 5 - 2, k 2 - 1, l 2 , l 3 , / 4 ; q) + q = y(kx -l,k 2 - 2, l 2 , l 3 , / 4 ; q) . 



(43) 



But, as explained above, the recursion with fixed indices l 2 ,l 3 ,k is not really enough since in the recursion for 
X 2 k+x, the boundary conditions will change, so the values for l 2 ,l 3 ,U will change correspondingly. We therefore 
derive a set of relations for the functions y(kx, k 2 ; l 2 , l 3 , k; q), which allow us to change the values of the U. 



2. Identities for y(ki, k2\ h, h, h', q) 



To derive the necessary identities, we start from the explicit definition of y(kx,k 2 ;l 2 ,l 3 ,k;q) in Eq. (40). The 
variables l 2 and /3 (both taking the values l 2 ,l 3 = 1,2) only appear in the g-binomials. The same is true for I4 if it 
takes the values 1 4 = 1,4. We can therefore relate the functions y(kx,k 2 ) for the two different values of l 2 (keeping 
l 3 ,k fixed), by shifting the values of kx,k 2 , and similarly for l 3 (with l 2 ,k fixed). To relate the functions with the 
values I4 — 2, 3, we need to swap the values of kx and k 2 and shift them, where the shifts depend on the values of 
l 2l l 3 . In particular, we find (suppressing the variable q, as we will frequently do as well below) 

y(kx, k 2 ; 1, l 3 ,l 4 ) = y{kx,k 2 + 1; 2, l 3 ,l 4 ) 
y(kx, fc 2 ; h, 1, h) = y{kx + l,k 2 - 1; l 2 , 2, Z 4 ) 
y(kx, k 2 ; l 2 ,l 3 , 1) = y(kx,k 2 ;l 2 , l 3 , 4) 

y(fci,fc 2 ;M,2) =y(k 2 ,kx; 1,1, 3) (44) 
y(kx,k 2 ; 1, 2, 2) = y(k 2 +2,kx~ 2; 1, 2, 3) 
y(kx,k 2 ; 2, 1, 2) = y(k 2 -l,kx + 1; 2, 1, 3) 
y(kx,k 2 ; 2, 2, 2) = y(k 2 + l,kx- 1; 2, 2, 3) . 



There are no further relations needed with I4 since U = a (see Sec. VA) and the reflection a 
change possible we can make in the recursion for X^+x- 



a is the only 



In view of the relations (44), there are two independent functions, say, 

y(kx,k 2 ; 1, 1, 1) and y(kx,k 2 ; 1, 1,2), 

which in particular satisfy the following identities 

y(kx,k 2 ; 1,1,1) = y(k 2 , kx; 1, 1, 1) 
y(kx,k 2 ; 1,1,2) = y(k 2 , kx; 1, 1, 3). 

In the following, we show that in fact y(kx, k 2 \ 1, 1, 2) = y(kx, k 2 ; 1, 1, 3) identically, without swapping the arguments 
ki,k 2 . From the point of view of the definition, Eq. (40), this is a rather nontrivial equation, because one can not 
simply relate terms in the sum of y(kx, k 2 ; 1, 1, 2) to terms in the sum of y(fci, k 2 ; 1, 1, 3); all the products of g-binomials 



get mixed. This is in contrast to the identities in Eq. (44), which could be obtained by trivial relabelings. Below, we 



derive similar equations for the other possible values of l 2 and l 3 . 

The recursion for y(kx, k 2 ; l 2l l 3 , 14) does not involve the variable l 4 . Using the recursion, it is easy to see that the 
initial conditions y(0, 0; l 2 , l 3 , 14), y(0, 1; l 2 ,l 3 , 14), y(l, 0; Z 2 , l 3l I4) and y(l, 1; l 2 , l 3 , 14) specify the values of the function 
y(kx, k 2 ; l 2 , l 3 , 14) uniquely. In fact, one can see that 



12 





(fa, fa) = 


(1,1) 


(1,2) 


(2,1) 


(2,2) 


y(0,0;fa,fa,2) 


= 2/(0,0, fa, fa, 3) = 


1 


1 /2 


1 


1 


2/(l,0;fa,fa,2) 


= 2/(1,0, fa, fa, 3) = 


1 + 9 


1 + 9 


g l/2 


1 


2/(0,1; fa, fa, 2) 


= 1/(0,1, fa, Is, 3) = 


1 + 9 


1 + 9 


1 


gl/2 


2/(M;fa,fa,2) 


= 2/(1, 1, fa, fa, 3) = 


2 gl/2 + ^3/2 


1 + 9 + 9 2 


1 + 9 


1 + 9 • 



This shows that the the functions y(ki, fc 2 ; fa, fa, 2) and y{k\, fc 2 ; fa, fa, 3) are in fact identical. 

Similar relations for Z4 — 1,4 are trivially true by changing the symmetric summation variables. So we see that, for 
r = 5, 



2/(fei,fc 2 ;fa,fa,fa) = y{kuk 2 ;l 2 ,h,r - fa). 



(45) 



These relations are analogous to the identity (26) for the height probability X 2 k+i(a; fa c, d, e). 
Using these with the last four relations in (44 1 gives the following, nontrivial relations 



2; 1,2, 2) 



y(h,k 2 ; 1,1, 2) = v (*a,A i; 1,1,2) 
y(k 1 ,k 2 ;l,2,2)=y(k 2 + 2,k 1 
y(h,k 2 ; 2, 1, 2) = y(k 2 - 1, fa + 1; 2, 1, 2 
2/(fa, k 2 ; 2, 2, 2) = y(k 2 + 1, fa - 1; 2, 2, 2 



(46) 



With help of Eq. (45), these equations also hold for Z4 = 3. In addition, by making use of the relations in Eq. (44) 
and shifting l 2 and fa when necessary, we find that they also hold for — and hence for Z4 = 4. Thus, we have 



y(h,k 2 ; 1, 1, fa) = y(k 2 , fa; 1, 1, fa) 
y(ki, k 2 ; 1, 2, h) = y(k 2 + 2, fa - 2; 1, 2, Z 4 ) 
2/(fa, fa; 2, 1, fa) = y(fc 2 - 1, fa + 1; 2, 1, l 4 ) 
2/(fa, fa; 2, 2, fc) = y(fc 2 + 1, h - 1; 2, 2, Z 4 ) 



(47) 



It is important to note that the functional form of the identities we derived does not depend on the value of fa. This 
will be very useful in the following, because the form of the recursion relations for y(kx, k 2 , fa, fa, fa) does not depend on 
fa either. In establishing the connection between X 2 ^ + i(a; 6, c, d, e) and y{k\, k 2 , fa, fa, fa), the cases which only differ 
in the values of fa (or a, which is the corresponding variable in the functions X 2 u+i) can be dealt with simultaneously. 



B. The identifications between X2k+i and y(k) 

We are now ready to state the relations between the functions X 2 k+i(a; b, c, d, e; q) and y(k±, k 2 ; fa, Z3, fa; q). Note 
that due to the properties discussed in sec. |II E[ these identifications have slightly different but equivalent form for 
the states with u < as compared to [5]. Using these identifications, we can show that the recursions for the functions 
y{ki, k 2 ) at special values of the arguments imply those of the functions X 2 k+i and that the initial conditions agree. 
Thus we are able to give the functions X 2 k+i in a closed form. 

We have two independent ground state patterns Gf and ■ For the ground states the identifications are 

X 2fe+1 (l;2,l,2,3 ;(7 ) - y(fc,fc;l,l,l;9) 

X 2fe+ i(3;2,l,2,3; g ) = y{k, k; 1, 1, 3; q) 

X 2fe+1 (2;3,2,3,4 ;( ?) = y(k, k; 1, 2, 2; q) 

X 2fe+1 (4;3,2,3,4; g ) = y(k, k; 1,2, 4; q). (48) 



and for the ground states G~{ 



+ 



X 2fe+1 (2;l,2,3,2;g) = y(k, k; 2, 1, 2; q) 
X 2fc+ i(4;l,2,3,2;g) = y{k, k; 2, 1, 4; q) 



k+1 



X 2fe+1 (l;2,3,4,3 ;(? ) = y(k,k;2,2,l;q) + q^y{k-l,k-l;l ) lA;q) 



k+1 



X 2fe+1 (3;2,3,4,3; g ) = y{k, k; 2, 2, 3; q) + q^y{k - 1, k - 1; 1, 1, 3; q). (49) 



These are all the (r — l)(r — 3) = 8 independent functions X 2 k+i(a; b, c, d, e; q) for r = 5. 
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From these identifications, all the ground state patterns specified in Ref. [5] are obtained using the properties of 
X 2k +\S and y's. In particular, the identities we derived here show that the identification of ground states 
is consistent with the somewhat different identification of ground states in terms of the function y(kx, k 2 ; hi h> h) m 
Ref. [5] for r = 5, because we always have that k\ = k 2 = k in the expressions for the height probabilities. 

We are now ready to show that the identifications made above are correct. We will first show that assuming that 
the identification is correct for X2k-i, this implies that the identification is also correct for X 2k+ \. To do this, we 
also make use of the recursion relations for X 2k+ i(a; 6, c, d, e; q) and the functions y{k\, k 2 ; h, h, 14', q), as well as the 
various relations between the latter. We will complete the proof by showing that the initial conditions also agree. 



1. Recursions for G 2 

To show that the identifications are indeed as given above, we frequently make use of the various relations between 
the functions y(ki, k 2 ; hi h', q) as given in section III A 2 For clarity, we repeat these relations when we use them, 
with one exception. We frequently use Eq. (45) to swap I4 O 5 — Z4, without mentioning this explicitly. 

We start by considering the cases G 2 , namely the recursion for X2fc+i(l; 2, 1, 2, 3) implies 

X 2k+1 (l; 2, 1, 2, 3) = q^ X 2k _ 1 (l; 2, 1, 2, 1) + X 2k _ 1 (l; 2, 3, 2, 1) + X 2k _ 1 {l; 4, 3, 2, 1) (50) 
or writing the RHS in terms of the functions y, 

X 2k+1 (l; 2, 1, 2, 3) - q^q k l 2 y{k - 1, k - 1; 1, 1, 1) + y(k —l,k— 1; 1, 2, 4) + q k ' 2 y{k - 1, k - 1; 2, 1, 4) . 
Now, y(k - 1, k - 1; 1, 2, 4) = y(k - 2, k; 1, 1, 1) and y(k - 1, k - 1; 2, 1, 4) = y(k - 1, k - 2; 1, 1, 1), so in total 
X 2k+1 (l; 2, 1, 2, 3) = q^q k ' 2 y{k - 1, k - 1; 1, 1, 1) + y(k - 2, fc; 1, 1, 1) + q k ' 2 y{k - 1, k - 2; 1, 1, 1) . 
Using the recursion y(k + 1, k — 2; 1, 1, 1) = q h i 1 y(k — 1, k — 1; 1, 1, 1) + y(k — 1, fc — 2; 1, 1, 1), we get 
X 2fe+ i(l; 2, 1, 2, 3) = q fc / 2 y(fc + 1, k - 2; 1, 1, 1) + y(k - 2, fc; 1, 1, 1) = y(k, k; 1, 1, 1) , 

where we first used the relation y{k\, k 2 \ 1, 1, 1) = y(k 2 , fei; 1, 1, 1), followed by the recursion relation for y(fc, k; 1, 1, 1). 
We have thus shown that indeed X2fe+i(l; 2, 1, 2, 3) = y(k, k; 1, 1, 1), based on the identification for 2k — 1 and the 
recursion relations. The case X 2k+1 (3; 2, 1, 2, 3) = y(fc, fc; 1,1,3) follows automatically, because all the relations we 
used are independent of the actual value for Z4 or a. 

The only G 2 case left to consider is X 2k+1 {2; 3, 2, 3, 4). The recursion for X 2k+1 {2\ 3, 2, 3, 4) is 



X 2k+1 (2; 3, 2, 3, 4) = q~ X 2k ^(2; 1, 2, 3, 2) + q~ X 2k ^{2: 3, 2, 3, 2) + X 2fe _ 1 (2; 3, 4, 3, 2) . 
Writing the RHS in terms of y's gives 

X 2fc+1 (2; 3, 2, 3, 4) = q^y(fc - 1, fc - 1; 2, 1, 2) + q^q k / 2 y(k - 1, fc - 1; 1, 2, 2) 

+y(fc- l,fc- 1;1,1,3) . 

Transforming everything to equal indices l 2l h,l4, we get 

X 2k+l {2; 3, 2, 3, 4) = q^ 1 (y{k, k - 3; 1, 2, 2) + q k ' 2 y(k - 1, fc - 1; 1, 2, 2)) + y(k, fc - 2; 1, 2, 2) 



(51) 



Now, eq. ( 46 ) implies 



Using this and the recursion gives 



y(h,k 2 ; 1, 2, 2) - y(k 2 + 2, fc x - 2; 1, 2, 2) . 



X 2fc+1 (2; 3, 2, 3, 4) - <? — y(fc + 1, fc - 2; 1, 2, 2) + y(fc, fc - 2; 1, 2, 2) 

which again just the recursion for y(k, fc; 1, 2, 2) = X2jt+i(2; 3, 2, 3, 4). The recursion for X2fe+i(4; 3, 2, 3, 4) follows 
similarly independent of I4 = a, using the identity |47|). 
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2. Recursions for G± 

For the ground states Gf, the first recursion for X 2 fc + i(2; 1, 2, 3, 2) leads to 

X 2fc+1 (2;l,2,3,2)=X 2fe _ 1 (2;3,2,l,2)+X 2fc _ 1 (2;l,2,l,2) . (52) 

In terms of the function y, the RHS is 

X 2fe+1 (2; 1, 2, 3, 2) = y(fc - 1, fc - 1; 2, 2, 3) + g fc / 2 y(fc — 2, fc — 2; 1, 1, 3) 

+flf fc y(fc-l,fe-l;2,2,2) 

or 

X 2fc+1 (2; 1, 2, 3, 2) = y(fc - 2, fe; 2, 1, 2) + <? fe / 2 y(fc - 2, fe - 1; 2, 1, 2) + q k y(k - 2, fc; 2, 1, 2) . 



Now ( 46 ) implies that 

y{k u fc 2 ; 2, 1, 2) = y(fc 2 - 1, ki + 1; 2, 1, 2) 

and we get, using the recursion once, 

X 2fe+1 (2; 1, 2, 3, 2) - y(k - 2, fe; 2, 1, 2) + q k ' 2 y{k - 2, fc + 1; 2, 1, 2) (53) 

which is again just the basic recursion. The recursion for X 2 fc +1 (4; 1, 2, 3, 2; g) follows similarly. 
The recursion for X 2 fc+i(l; 2, 3, 4, 3; q) leads to 

X 2k+1 (l; 2, 3, 4, 3) = g^X^l; 2, 1, 2, 3) + X 2fe _ 1 (l; 2, 3, 2, 3) + X 2fe _i(l; 4, 3, 2, 3). (54) 

Writing this in terms of y's is 

X 2fe+1 (l; 2, 3, 4, 3) = g^y(fc - 1, fe - 1; 1, 1, 1) + q k ' 2 y{k - 1, fe - 1; 1, 2, 1) + y(fe - 1, fc - 1; 2, 1, 4). 

This is simply 

X 2fe+1 (l; 2, 3, 4, 3) = q^y{k - 1, fc - 1; 1, 1, 1) + <z fc / 2 y(fc - 1, fc; 2, 2, 1) + y(k, k - 2; 2, 2, 1). 



Now eq. ( 47 1 implies 

y(fc x , fc 2 ; 2, 2, 1) = y(fc 2 + 1, fci - 1; 2, 2, 1), 

using this and the recursion gives back y(k, fc; 2, 2, 1) + q h ^y{k — 1, fc — 1; 1, 1, 1) = X 2 fc +1 (l; 2, 3, 4, 3). The recursion 
X 2 fc + i(3; 2, 3, 4, 3) is identical. 

Since Gf and give all the independent states and recursions for X^k+x, we have shown that the recursion for 
y(fci, fc 2 ; h, h, h) implies the recursion for X 2 fc_|_i(a; b, c, d, e). 

3. Initial conditions 

To establish the equality of the functions X 2 j, +1 (a; b, c, c?, e) and the functions y(fc, fe, / 2 , Z 3 , Z4), we still have to verify 
the initial conditions. We first show that knowing the functions y(fci, fc 2 ; Z 2 , 13, 14) for fci, fc 2 = 0, 1 for all values of the 
li fixes the functions completely, by means of the recursion relations. We include the argument in detail here, because 
we will need it in the discussion of the case r = 6. 

It follows from the recursion relations that if we know the functions y for all (fci, fe 2 ) with k\ + fc 2 < n, where n > 2, 
we can construct all the functions with k\ + fc 2 = n + 1. In the following, we will suppress the dependence on Z 2 , Z3, 14. 
Namely, the function y(i, n — i + 1), where i < n/2, can be obtained from y(i, n — i — 1) and y(i + 1, n — i — 1), which 
we both know by assumption. Similarly, we can obtain y(n — i + from y(n — i — 1, i) and y(n — i — l,i + 1). 

From the knowledge of the functions y(0, 0), y(0, 1), y(l, 0) and 1), we first obtain y(0, 2) and y(2, 0) from the 
recursion. The argument above shows that we now can obtain all the functions. 

It is now a simple matter to check that the initial conditions for the functions X 2 fc_)_i(a; b, c, d, e), namely 
Xi(a; 6, c, d, e) and ^3(0; b, c, d, e), indeed correspond to y(0, 0; Z 2 , 13, 1 4) and y(l, 1; i 2 , 13, 

Below, we give these initial conditions for the independent set of functions Xzk+x, namely those corresponding to 
the ground state patterns for u > 0. For the ground state patterns of type G 2 , we have 
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(a, b) 


(1,2) 


(2,3) 


(3,2) 


(4,3) 


Xi(a; 6, 6 — 
X 3 (a; 6, 6 - 


1,6,6+1) 
1,6,6 + 1) 


gl/2 
l + q 2 


1 + q + q 2 


1 

2q l/2 + ^3/2 


1 

g l/2 + g 3/2 



The results for the ground state patterns of type read 



(a, 6) 


(1,2) 


(2,1) 


(3,2) 


(4,1) 


Xx (a; 6, b + 1,6 + 2, 6 + 1) 
X 3 (a;6,6+l,6 + 2,6 + l) 


qi/2 

(Z l/2 +g 3/2 


1 

1 + 9 


1 

l + 2g 




q i/2 



This shows that the initial conditions are identical and completes the proof. 



IV. EXPLICIT EXPRESSIONS FOR r = 6 

We continue by considering the case r — 6. The form of the functions y( A,B ' C \kx, fc 2 , ^3; ^2, '3, ^4; <?) will be motivated 
in section [Vj where we give a conjectural form for the functions y for general r. 



A. The function y (A ' B ' c) (ki, k 2 , k 3 ; l 2 , h, h\ q) for r = 6 
For r — 6, the 'finitized CFT characters' related to the functions X 2 fe+i(a; 6, c, d, e; (?) take the form 



t A ' B < c Hkx,k2,k 3 ;l 2 ,l 3 ,h;q) = 



-5(mi+m2+m|— (mi+m3)m2— mii5( 4l 4— OT2<5! 4 ,3— TO35( 4 ,2) 



(mi,m 2) m 3 )6(2Z> ) 3 + (A,B,C) mod 2 



fci+m 2 +<5i, ,1 +<5 



2 

mi 



fc3+™l+™3+<5; 3 ,2+t>"l 4 ,3 

2 

m 2 



fc 2 +m2+l?i 2 ,l+l?i3,3+^4, 

2 

m 3 



(55) 



where the summations are restricted such that the rrii have the same parities as (A, B,C), which now have to be 
explicitly specified in contrast to the case r = 5. We will only consider these functions for (i 2 + l 3 + I a) mod 2=1. 
The properties of this function derived below are valid under this condition. 

We start by noting that the condition (7 2 + h + ^4) mod 2 = 1 implies that for the function y r= e to be non-zero, 
one needs that fci = fc 2 mod 2. This follows from the requirement that the arguments of the g-bmomials have to be 
integers. Thus, one requires that both B = (fci +<^ 3 , 1+^4,4) mod 2 and B = (fc 2 + <5z 2 ,i+<5z 3 ,3 + <5z 4l2 ) mod 2. Inspection 
shows that to satisfy both equations, one needs fci = fc 2 mod 2. In addition, one needs that A + C — (k 3 + <^ 3 ,2 +^4,3) 
mod 2, so there are two, a priori independent, functions for every £2,^3, '4, given by the two different choices for A 
and C. 

Again, the height probabilities X2k+i are related to these more general functions with kx = fc 2 , k 3 = 0. The 
variable £3 is introduced in y( A,B,c \kx, k 2 , ^3) to obtain a recursion that closes. When we state the relation between 
the functions X 2 ^. +1 (a; 6, c, d, e; q) and y( A,B ' C \k, k, 0; Z 2 , l 3 , 14; q) below, we will specify the required values of (A, B, C) 
explicitly. 



1. Recursion for y^ ABC \ki, fe, £3; h, h, h', q) 



The recursion for this function, following from (|39| is, with I = l2,h,h to shorten the notation, 

y 



( A ' B > C) (kx,k 2 , fc 3 ; l; q) = q f^y( A +^ B ^ (kx - 2, fc 2 , fc 3 + 1; £■ q) + y^ c \kx - 2, fc 2 , fc 3 ; 4 <?) 



1 WJ ' ^(fci, fc 2 - 2, fc 3 + l;4 9) + y (A ' B ' c) (kx, k 2 -2,k 3 -J;q) (56) 

(A ' s ' C) (fci,fc2,fc 3 - 2;£; q). 



= q 2 y 

= q 3+Sl t 2 " y(A,B+X,C) ^ +1>k2 + 1; ^ _ 2 ; q) + y 



As was the case for r — 5, we need to establish several identities for the functions y^ A ' B,c \kx, k 2 , k 3 ; i 2 , '3, h', 9), 
which allow us to change the values of I2, l 3 , 14, in order to make connection with the recursion for X2k+i, as given in 
equation ( 25 ) , which does not keep the U constant throughout the recursion. 
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2. Identities for y^ A ' B ' c \k\, ki, k%; h, h, h; q) 

The relations we need in order to prove the equivalence between the X2k+i and y^ A,B,c '{k\, k 2 , k 3 ) can be derived 
in the same way as for r = 5. However, we need to keep track of the parities (A, B, C), complicating matters slightly. 
As explained above, we only consider the cases obeying (7 2 + l 3 + h) mod 2 = 1. 

Similar to the situation for r = 5, the relations can be grouped into several classes. We first deal with the relations 
which can be obtained from the definition ( 55 1 by shifting the values of the fc, , and if necessary changing the summation 
variables m\ f-> m 3 . All these relations are therefore trivial in nature. 

The first class of relations relates the functions which only differ in the values of l 2 and l 3 . For l 4 odd, we have 
(dropping the q dependence) 

(h, k 2 , k 3 ; 1, 1, U) = y( A ' B ^ (h + l,k 2 + 1, k 3 - 1; 2, 2, l 4 ) = y^ c \h + l,k 2 - 1, h; 1, 3, h) . (57) 

In the case I4 even, we find 

y^ c \h, k 2 , k 3 ; 1, 2, h) = y^'^ih -l,k 2 + 1, k 3 + 1; 2, 1, l 4 ) = y^ B ^ c \k u k 2 , h + 1; 2, 3, Z 4 ) • (58) 

We continue by relating the functions y with l 4 and r — l 4 , by changing the summation variables mi ■<-> 7713, which 
also swaps the values of ki and k 2 . For Z 4 odd, this gives 

y (A ' B ^(ki, k 2 ,k 3 , 1, 1, l 4 ) = y {c - B < A \k 2 , k u k 3 , 1, l,r - k) 
y^ B ' c \k x , k 2 ,k 3 ,2, 2, l 4 ) = y ic ' B ' A) (k 2 , k u k 3 , 2, 2,r - J 4 ) 

^■^(fci, fc 2 , /c 3j 1, 3, 1 4 ) = y^ B ' A \k 2 + 2, fci - 2, fc 3 , 1, 3,r - / 4 ) • (59) 

For / 4 even, we find 

y{A,B,c) ^ fcaj fc3) 1; 2j j 4 j = y (C,fl,A) (fca + 1: fcl _ 1; ^ l2r _ h) 
y (A ' B ^(ki, k 2 ,k 3 ,2, 1, h) = y {c - B > A) (k 2 - 1, kx + 1, k 3 , 2,1, r- U) 

y( A ' B < c \k u k 2 ,k 3 ,2, 3, h) = y^ B ' A \k 2 + 1, k x - 1, fc 3 , 2, 3,r - l 4 ) . (60) 

We now relate the functions y^ A ' B,c >{k\, k 2 , k 3 ; l 2 , l 3 , £ 4 ) with l 4 and r — l 4 , but without changing the summation 
variables mi -f-> 7713. This can be done trivially in the case I4 — 1,5, but in the other cases, the relations are nontrivial, 
because they completely scramble the contributions from the products of the binomials, and hence can not be obtained 
by reshuffling the terms in the sums in equation (55). Instead, these relations arc obtained by using that the recursion 
relations for y are independent of I4, and checking the initial conditions. The initial conditions sometimes give rise to 
constraints for the ki. We always assume that ki,k 2 > 0, and give the constraint on £3 explicitly. 

y( A > B 'V (ki,k 2 , k 3 , l 2 , l 3 ,l) = y( A - B >V (ki,k 2 , k 3 , l 2 ,l 3 , 5) 

y^ A ' B ' c \h,k 2) k 3 , l 2 , l 3 ,3) = y {A+1 > B > c+1) (ki,k 2 , k 3 , l 2 ,l 3 , 3) for k 3 + S h - 2 > 

y^ B ^(ki,k 2 ,k 3 ,l 2 ,l 3 ,2)=y^ A+1 - B +^ c+1 \ki,k 2 ,k 3 ,l 2 ,l 3 ,A) for k 3 + 5 l3 , 2 > . (61) 

Finally, we combine these identities with the preceding ones, to obtain expressions relating the functions y with the 
same values of l 2 ,l 3 ,h, but with the values of k\ and k 2 swapped. In particular, for l 4 odd we have 

y^^Hk^k^h, 1, 1, h) = y( c + 5 '^> B ' A + s U,*)(k 2 , ki,k 3 , 1, 1, Z 4 ) for k 3 > 

y^ B ' c Hki, k 2 , k 3 , 2, 2, l 4 ) = y<. c +>U.* B . A +*U.»)fa, h,k 3 , 2, 2, l 4 ) for k 3 > -1 

y^ A < B ' c \ki,k 2 , k 3 , 1, 3, h) = y( c + s ^> B > A + 5 H,*)(k 2 + 2,h- 2, k 3 , 1, 3, Z 4 ) for k 3 > . (62) 

For l 4 even, we finally obtain 

y {A - B ^ (ki,k 2 , k 3 , 1, 2, l A ) = y (C+hB+i,A+i) (A . 2 + ljkl _ 1; fc3j 1; 2] u) for k 3 > -1 

y^ B ^ (ki ,k 2 ,k 3 ,2,l,l 4 )= y tc+i.B+i,A+i) (fe 2 - 1,^ + 1,^3,2,1,^4) for fc 3 > 

y( A ' B ^ (ki ,k 2 ,k 3 ,2, 3, U) = y (C+hB+i,A+i) (^ + 1,^-1,^,2,3, h) for fc 3 > . (63) 

This exhausts the relations that we need to prove the equivalence between the functions X 2 k+i and y( A ' B > c } (fci, k 2 , k 3 ). 
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B. The identifications between X2k+i and y^ A ' B,c \ki, fe, fe) 



As was the case for r = 5, it suffices to give the identifications for an independent set of X2k+i- We again specify 
the cases corresponding to the ground state patterns G\~ and G^ . We henceforth drop commas from the boundary 



conditions, b, c,d,e — > bede, the parities of the summation variables A, B, C 
to lighten the notation. 

For the patterns in G^~, the identifications are 



ABC and the labels ^2,^3,^4 



^2fc+l(l 
^2fc+l(3 

^2fc+l(l 
-^2fc+l(3 

X 2 k+i(5 

In the case of patterns of type , they read 



2123) 
2123) 
2123) 
3234) 
3234) 
4345) 
4345) 
4345) 



= y^)(fc,fc,0;lll) 
= ^ c )(fc,fc,0;113) 

= y<^>(fc,M;115) 
= y^)(fc,fe,0;122) 



OL4C) 



(fc,fc,0;124) 
y ^)(fc,fc,0;131) 



y 



(ABB) 



(fc,fe,0;133) 
y^ AA )(fe,A;,0;135). 



(64) 





.1 (2; 1232) 


= y^ AAA \k,k,Q;2l2) 












,i(4; 1232) 


= / ABA \k,k,0-,m 










-^2fc- 


-i(l;2343) 


= y( ASB )(fc,fc,0;221)- 


fe+1 

\-q 2 


y^ BBB \k- 


l,k- 


1,0; 111) 


X2k- 


_!(3;2343) 


= y (ASA) (M,0;223)- 


fe+i 
Vq 2 


y( BBA \k~ 


l,k- 


1,0; 113) 


X2k- 


_i(5; 2343) 


= y (AAC \k,k,Q- 22b) - 


k+i 
Vq 2 


yt CAC \k- 


l,k- 


1,0; 115) 


X2k- 


-1 (2; 3454) 


= y (ASA) (M,0;232)- 


fe+i 
Vq 2 


y (BAA \k- 


l,k- 


1,0; 122) + q^y 


X2k- 


-1 (4; 3454) 


= y (AAA) (fc,fc,0;234)- 


fc+i 
-q — 


y( BBA \k- 


l,k- 


1,0; 124) + q^y 



( BBB \k- l,k- 1,0; 212) 
(BAB \k- l,k- 1,0; 214). (65) 

We still need to specify the parities (A, B, C) for the functions y, in order to completely determine the identification. 
These read as follows 



A = {k + l 3 + S k>5 + ^3,2) mod 2 
B = (k + l 3 + S hA + ^3,3) mod 2 
C = (k + l 3 + <5 U)3 + <y Ui5 ) mod 2 



(66) 



It follows that all the parities are reversed whether considering k even or odd. The notation for the parities in the 
equations (65) and (64) requires some explanation. Clearly we need at most two different labels for the parities. The 
parity A of the first function y( A ' B,c ^ (k, k, 0; 1%, I3, 14) for size k is used as a reference, and is always denoted by A. 



This parity is obtained from equation (66). If the parity for B, also given by equation (66), happens to be the same 
as A, this is also denoted A. Otherwise, it is denoted as B. Finally, for the parity C, the notation is such that it is 
denoted by A if A = B = C, denoted by B if A ^ B = C, and denoted by C if A = B ^ C. When the value of k 



is lowered, as occurs for the Gf patterns, see equation (65), the notation of the parities is with respect to those of 
size k. This notation is convenient to keep track of the parities, when we prove (in the next section) the connection 
between X2k+i and the y^ A ' B - c \k, k, 0), as given in (65) and (64), and write the parities with respect to the function 
for size k as described above. 

The identifications given above deal with all the (r — l)(r — 3) = 15 independent configurations for the functions 
X2k+i for r = 6. The relations for the functions X 2 k+i for the ground state patterns related to u < 0, and the non 
ground state patterns, can be obtained by making use of the relations given in the equations (32), (33) and (31). 



1. Recursion for states G~t 

Here we show how to deal with a representative state in Gf , the rest of the recursions are collected in appendix [a] 
The strategy is the same as for r = 5, namely we assume that the identification is correct for 2k — 1, and show that 
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this implies the identification for 2k + 1, by using the recursions for X 2k+ \ and y^ A,B ' C \k\^k 2 ^k 3 ) as well as the 
various relations between the functions y. The recursion for X 2 fc + i(3; 2343) is 

X 2fe+1 (3;2343)=5** i X 2fc _ 1 (3;2123) + X 2fe _ 1 (3;2323)+J5!: 2fe _ 1 (3;4323) (67) 
In terms of y^ ABC \k 1 ,k 2l k 3 ;l 2> h, h) the RHS is 

X 2k+l (3; 2343) = q^y^ BBA \k - 1, k - 1, 0; 113) + qiy^ BAA \k - 1, k - 1,0; 133) 
+y {BAB \k - 1, k - 1,0; 223) + q^y {AAB \k - 2, k - 2, 0; 113). 

We use 

y( ABC )(A :i ,fc 2 ,fc 3; 133) - 2 / ABC )(fe 1 ,fc 2 + 2,/c 3 -l;223), 
y( ABC )(A :i ,fc 2 ,fc3;113) = y (ABC) (fci + l,*c 2 + l,* 3 -l; 223), 

except for the first term which is included in X 2k+ i(3; 2343), to get 

X 2fe+ i(3; 2343) = q"^ 1 y {B ' B ' A) (k - 1, k - 1, 0; 113) + y {BAB) {k - 1, k - 1, 0; 223) 
+g V SAA >(fc - 1, * + 1, -1; 223) + q V AAS )(fc - 1, A: - 1, -1; 223). 
Where the term the second term on the RHS is 

y(BAB) ( fc _ i : A: - 1,0; 223) = y {BAB) {k - 1, k - 1, -2; 223) + (fc, fc, -2; 223), 

so 

X 2k+1 (3; 2343) = q^y^ 3 ^ (fc — l,fc — 1,0; 113) + y (Bj4S) (fc - 1, fc - 1, -2; 223) + (fc, k, -2; 223) 

+ g ly(- B ^)(fe - 1, k + 1,-1; 223) + qiy {AAB \k — 1, fc — 1, -1; 223). 

Combining the second and last terms gives 

X 2k+l (3; 2343) = q^y^ BBA ) {k - 1, k - 1, 0; 113) + j/ SAB ) (fc + 1, fc - 1, -2; 223) + j/ 53 ^ (fc, fc, -2; 223) 

+qiy( BAA \k-l,k + 1,-1; 223) 



and further, using the properties in (62) and (59), 

y( ABC Hk 1 ,k 2 ,k 3 ;223) = ^+ 1BC+1 )(fc 1 , fc 2 , fc 3 ; 223) = y( CSA > (fc 2 , fc l7 fc 3 ; 223), k 3 > -1, 

in the last term gives 

X 2k+1 (3; 2343) = q^y( BBA ^ (k - 1, k - 1, 0; 113) + 2/ (Bj4S) (fc + 1, fc - 1, -2; 223) + ?/ SBB) (fc, fc, -2; 223) 

+ q ^ y ( BAA )(k + l,k~ 1,-1; 223). 

This is just 

X 2fe+1 (3; 2343) = q^y^ 3 ^ (fc - 1, k - I, 0; 113) + y< BXfl > (fc + 1, fc + 1, -2; 223) + y( BBB ^ (fc, fc, -2; 223) 

= q h ^y^ BBA \k - 1, fc - 1, 0; 113) + y^ BBB Hk, fc, 0; 223) 
= q ^ y i BBA ) (fc - 1, fc - 1, 0; 113) + y {ABA) (k, fc, 0; 223), 

as desired. The recursions for the same patterns but different a — I4, i.e. X 2 fc + i(l; 2343) and X 2k+ i(5; 2343), are 
similar and omitted. The remaining recursions for patterns in are collected in appendix |a} 
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2. Recursion for states G 2 
Here we show an example recursion for a state in G 2 ■ We start with 

X 2k+1 {3; 2123) = q*Px 2 k-i(3', 2121) + X 2k _ 1 (3; 2321) + X 2 fc-i(3; 4321). (68) 
In terms of the functions y^ ABC \k\, fc 2 , k^; l 2 , Z 3 , I4), the RHS is 

X 2k+1 {3; 2123) = g^g V CCA) (k - 1, k - 1, 0; 113) + yl CAA > {k - 1, k - 1, 0; 133) 
+g t y (CVlC) ( fc _ fc _ i ; 0; 223) + q k y( AAC ) (k - 2, fc - 2, 0; 113). 



Using 



this is 



^ BC )(fc 1; fc 2 ,fc 3 ;223) = ^ sc )(fci ~ 1, ^ - 1, fc 3 + 1; 113), 
2/ (ABC) (fci,fc 2 ,fc 3 ;l33) = ^ sc Hfci-l,fc2 + i,fc 3 ;H3), 



X 2fc+1 (3; 2123) = q^qiy {CCA \k - 1, A - 1,0; 113) + j/ CAA >(fc - 2, fc, 0; 113) 
+qiy^ CAC \k — 2, fc — 2, 1; 113) + q k y( AAC \k - 2, fc - 2, 0; 113). 



We use 2/ ASC, )(fci, fc 2 , fc 3 ; 113) = y^ CBA \k 2 , fci, fc>3; 113) from (59) and combine it with the recursion, 

y(AAC) ( fc _ 2, - 2, 2; 113) = q l / 2 y( CCA ) (fc - 1, fc - 1, 0; 113) + 2/ AAC ) (fc - 2, fc - 2, 0; 113), 



and get 



Using 



we get 



X 2k+1 (3; 2123) = gV AAC )(fc - 2, fc - 2, 2; 113) + y^ CAA \k - 2, fc, 0; 113) 

+9 fc/2 y (CAC) (A ._ 2)fe _ 2)1;113) 

y {CAC \k, fc - 2, 1; 113) = j/ CAC )(fc - 2, k - 2, 1; 113) + qiy {AAC ^y{k - 2, fc - 2, 2; 113), 
X 2fc+1 (3; 2123) = <zV CAC) (fc» - 2, 1; 113) + y (CAA) (fc - 2, fc, 0; 113) = y (AAC \k, fc, 0; 113), 



where we used y( CAA \k - 2, fc,0; 113) = j/ AAC )(fc - 2, fc,0; 113) (recall that the parity A f C) and the recursion for 
y. Thus, the identification X 2k+ i(3; 2123) = y^ AAC ^ (fc, fc, 0; 113) follows from the identifications for 2fc — 1 and the 
recursion for y, as we intended to show. 

The recursions for X 2 A:+i(l; 2123) and X 2k+ i(5; 2123) with different a — I4 are very similar and therefore omitted. 
The recursions for the other patterns in G 2 are given in appendix [X] 



3. Initial conditions 

We have now proven the recursions for the independent states in and G 2 are implied by the recursion for the 
function y r =6 and it remains to check that the initial conditions agree. 

We show that to specify the functions y r =e completely, it suffices to know them for all the values of ki £ {0, 1}. 
The others then follow by making use of the recursion relations. The argument is very similar to the one we gave for 
the case r = 5. 

We start by giving the recursion relations in symbolic form, 

(fci,fc 2 ,fc 3 ) -> (fci - 2, fc 2 ,fc 3 ) + (fci - 2, fc 2 ,fc 3 + 1) , 
and similarly for k\ -H- fc 2 , which occur symmetrically in the recursion. The value of fc 3 can be lowered by using 

(fci,fc 2 ,fc 3 ) -> (fci,fc 2 ,fc 3 -2) + (fc 1 + l,fc 2 + l,fc 3 -2) . 
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It is straightforward to show that from the initial values, we can obtain all functions y r =6 where the ki satisfy 
ki + fc 2 + fc 3 < 3. We start by constructing (^1,^2,^3) = (2,0,0), (0,2,0) and (0,0,2), by using the functions for 
the values (0,0,0), (0,0, 1) and (1, 1,0), which we know by assumption. We can also construct (3,0,0), (0,3,0) and 
(0,0,3), by using the values for (1,0,0), (1,0, 1), (0,1,0), (0, 1,1) and finally (0,0, 1), and (1,1, 1). The functions for 
these values also give us the functions at the values (2, 1, 0) and (1, 2, 0). 

Above, we constructed y r= § for (0, 0, 2), which together with (0, 0, 1) gives us y r= § for (2, 0, 1) and (0, 2, 1). Finally, 
we construct y r= § for the ki values (1, 0, 2) and (0, 1, 2) by using (1, 0, 0), (0, 1, 0) and the newly constructed (1, 2, 0) 
and (2, 1, 0). Thus, we now know y r =e for ah ki with k\ + fc 2 + k 3 < 3. 

We continue by showing that knowing the function y r =e for k\ + k 2 + ^3 < n allows us to construct the function 
for all ki values satisfying fci + k 2 + k% = n + 1. We start by using the first recursion, or its equivalent by swapping 
ki <-> k 2l to construct y, = 6 f° r the values (i, n — i + 1, 0), by using (i, n — i — 1, 0) and (i, n — i — 1, 1) or by using 
(i — 2, n — i + 1, 0) and (i — 2, n — i + 1, 1), whichever is applicable. We can similarly construct (z, n — i, 1) by using 
either (i, n — i — 2, 1) and (i, n — i — 2, 2) or (i — 2, n — i, 1) and (i — 2, n — i, 2). 

From now on, we can exclusively use the second recursion, to construct the remaining functions. We start by 
constructing (i,n — i — 1,2) from (i, n — i — 1,0) and (i + l,n — i,0) (which was just constructed above). We 
then continue to subsequently construct (i,n — i — j + 1, j), for j — 3, . . . , n + 1, from (i, n — i — j + l,j — 2) and 
(i + 1, n — i — j + 2,j — 2), which have been constructed earlier. This concludes the proof that we can construct 
the function y r =6 for all values of [k±, fc 2 , ^3) from the knowledge of the functions for the values ki € {0, 1} and the 
recursion relations. 

The initial conditions are now straightforward to check, namely the functions y^ ABC '{G, 0, 0; 1%, I3, 
y( ABC \l, 1, 0; I2, h, h) and X\(a; 6, c, d, e), ^"3(0; b, c, d, e) agree. Explicitly, the initial conditions read 



(0,6) 


(1,2) 


(1,4) 


(2,3) 


(3,2) 


(3,4) 


(4,3) 


(5,2) 


(5,4) 


Xi(a; b,b — 
X 3 (a;b,b~ 


1,6,6+1) 
1,6,6+1) 


l + q 2 




q 3/2 


qi/2 
1 + q + q 2 


1 

2q l/2 + g 3/2 


ql/2 
1 + q + q 2 


1 

2 5 l/2 + g 3/2 




q 


1 

q l/2 + q 3/2 



for ground state patterns of type G 2 , while for the ground state patterns of type we have 



(a, 6) 


(1,2) 


(2,1) 


(2,3) 


(3,2) 


(4,1) 


(4,3) 


(5,2) 


Xi(a;6,6+l,6 + 2,6+l) 
X 3 (a;6,6+l,6 + 2,6+l) 


q 1 ' 2 
q l/2 + q3 /2 


1 

1 + 9 


q 1 ' 2 
ql/2 + 2g 3/2 


1 

l + 2q 




q l/2 


1 

l + 2q 




ql/2 



This finally completes the proof of the co rrespo ndence between the functions X 2 fe + i(a; 6, c, d, e) and 
y( ABC \k, k, 0; I2J3, h) for r = 6, as given in section IVB 



V. CONJECTURAL EXPRESSION FOR THE LOCAL HEIGHT PROBABILITIES: GENERAL CASE 

In the previous two sections, we proved the equality between the explicit expressions (in terms of unitized CFT 
characters), and the functions X 2 k+\, in the cases r — 5 and r — 6. Now we present a general conjecture for X 2 k+i 
for any r > 5. We start this section by motivating the form of our conjecture. 

The form of the fermionic formula for general r closely resembles the fermionic formula for the height probabilities 
of the original Andrews-Baxter- Forrester model, which were first conjectured in [171 118) . and subsequently proven in 
several papers [HmHlS]. 

The structure of these fermionic characters is that of the 'Universal Chiral Partition Function' (UCPF), which can 
be interpreted as a sum of states of several species of fermions |13j . This formalism has been developed in the nineties 
by several groups, focussing on different aspects of the problem, such as the statistical mechanics models in which 
these UPCFs appear [53] . representation theory [TMTB] . or the connection with conformal field theory [331 155] , 

For an introduction and more references, see the note |13j . Here, we would like to mention in particular the paper 
|37j . in which the role of various types of particles which can appear in a UCPF is explained. By making use of the 
connection with the central charge, we can find the 'particle content' of the UCPF, which in turn dictates its detailed 
form. 

In the UCPFs which describe the local height probabilities on the one hand, and the characters of the CFT 
describing the critical behaviour on the other, two different types of particles appear, namely 'real' particles and 
so-called 'pseudo' particles, which in effect have zero-energy and do not propagate. Their presence does however effect 
the possible energies for the 'real' particles. In addition, the presence of 'pseudo' particles has an effect on the central 
charge associated with the UCPFs, which is why we bring up the notion of 'pseudo' particles here. 
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We start by considering a set of particles, satisfying Haldane exclusion statistics [35] , with the statistics parameters 
encoded in the matrix K, where the elements Kij denotes the mutual statistics between particles of type i and j. The 
matrix K plays a central role in the UCPF. The one-particle distribution functions A^ for an ideal gas of fractional 
statistics particles were derived in [5UH41) . 

where Zi = e~@ ei is the fugacity of the species i. Introducing A to t = Ili where all the Aj are evaluated at the 

same fugacity z, we can derive the central charge of the associated CFT by making use of the relation between the 
central charge and the specific heat. For the details, we refer to [37]. In short, the central charge associated with a 
system of exclusion statistics particles characterized by K (assuming that all particles are real for now) is obtained 
by calculating the specific heat. Specifically, one finds one has to solve the system of equations 



6 1 I ' 1 &)* • (70) 



The solution is then used to calculate the central charge 



c =^E L te), (71) 



where L(z) is Rogers' dilogarithm 



2 Jo \l-y v ) 

For pseudo-particles, which keep track of the internal structure of the real particles, we have ej = 0, so Z\ = 1, which 



changes the Aj via Eq. (69 1. To calculate the central charge in the presence of pseudo-particles, one first calculates the 
central charge associated with the full matrix K, giving C{ u \\. In addition, one calculates the central charge associated 
with the part of the matrix K which only contains the pseudo-particles, which we denote by c pseu <jo- The total central 
charge of the system is simply the difference of the two 

CCFT = Cf u n — Cp S0U do ■ (73) 

The statistics matrix K plays a central role in the UCPFs, which describe the partition functions of statistical 
mechanics models, as well as the (chiral) characters conformal field theories. For the various sectors in CFTs, both 
the matrix K and the type of the particles (real or pseudo) is the same. The forms of the functions y for r = 5 and 
r = 6 take the form of a UCPF. The bilinear forms appearing in the exponent of q, can be written as g2 mKm (see 



Eqns. (34 1,(55 1 for the cases r = 5 and r = 6). We use the discussion above to obtain the correct form of the matrix 
K, and the number of pseudo-particles necessary for the explicit expression for the functions X2k+i{a] b, c, d, e; q) for 
arbitrary r, but consider the original ABF model first. 

In the case of the original ABF model, we know that the critical behaviour of the model in regime III, i.e. when p->0 
is given in terms of the minimal models A4(r — l,r), which have a coset description s "( 2 ^*^M(2) r _ 3 w ^ n r > 4 |42j . 
The matrix K which appears in the bilinear form in the UCPFs for these minimal models was found to take the form 
K = |A r _3, where A r _3 is the Cartan matrix of su(r — 2), and is given via its elements as (A r — 3)1,3 = 2Si,j — 1, 
with i,j = 1, . . . , r — 3. To calculate the central charge from this matrix, wc need the additional information that only 
the particle associated with the first row and column is a real particle, all the other particles are pseudo-particles. 



We quote the results of the calculation here, and refer to j33j for the details. The solution to the equation ( 70 1 for 
the matrix K = |A. r _3 is given by 

sin(^) 2 



(1 — 2+j)n \ 2 



with j = 1, 2, . . . , r — 3. To obtain the central charge, dilogarithm identities were used, see [43], which resulted 
in the central charge associated with the matrix K = lA,._ 3 (assuming that all particles are real), namely c r = 

■^2 Sj=i = ( r ~ 3 - ) J r ~ 2 - ) . Thus, the central charge of the minimal models, for which only the first particle is real, 

is given by c mm = c r — c. r _i = 1 — ^rzfjz^ , which is indeed the central charge of the minimal model M(r — 1, r). We 
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also note the following. The sum of the central charges associated with a matrix and its inverse add up to the rank of 
the matrix (assuming that all particles are 'real'). This means that the central charge associated with K -1 = 2A~_^ 3 
is given by r — 3— ( r ~ 3 )J r ~ 2 ) = 2< " r ^ 3 ^ , which is the central charge associated with the Z r _ 2 -pa r afermion theory, which 
describes the critical behaviour of the critical p — model for u < 0, i.e. the critical behaviour associated with regime 
II of the ABF model. We explain this connection in more detail below. 

We now focus on the CFT description of the composite height model introduced in [5]. There, it was found that the 
critical behaviour in regime III was given in terms of a diagonal coset model, also based on su(2) r affine Lie algebras, 
in particular SM ( 2 )i x ^( 2 )i x sM(2) r _ 4 _ j n jjj^ explicit fermionic expressions for the local height probabilities, for finite 
system size, were found for the special case where r = 5. These local height probabilities took the form of UCPFs, 
based on the matrix K = ^A 2 , where both particles are real. The associated central charge is indeed C5 = |. 

The central charge of the coset sn(2)lXsn ^ lXstl(2) '- 4 is given by c = 1 + 1 + - = 2 - . B y using the 

analogy with the original ABF model, we can expect that the LHP of the composite model can be expressed in terms 
of UCPFs, based on the matrix K = ^A r _3, where two of the particles are real, while the others are pseudo-particles. 
Indeed, the central charge associated with such a UCPF is given by c r — c r _2 — 2 — r fZ_ 2 \ , which is the expected 
result. Below, we will give the conjectured form of the LHP, written in terms of fermionic characters, based on the 
bilinear form of K = |A r _3, where the first and last particles are real, while the other particles are pseudo-particles. 



A. Conjectured fermionic form of the local height probability 

In this subsection, we will introduce the conjectured form of the functions X r 2fe+i (<x; b, c, d, e; q) which appear in 
the expression for the local height probability, for arbitrary value of r, thereby generalizing the result given in [5] for 
r = 5, which we proved in the current paper, along with the case r = 6. Note that we added the subscript r to the 
notation X Tt 2k+i(fl\ b, c, d, e; q), in order to be completely explicit. 

We first introduce a set of functions y r (k; I2, I3, I4] q), which play a central role in the description of the functions 
X r .2fc+i(«; b, c, d, e; q). In particular, we write 



y r (k;l2,h,k;q) = 



r ■ 

mi>0 
=l,...r-3 



|(fe + m 2 + <5; 3 ,i + <5i 4 , r -2) 
mi 



r— 4 

n 



+ m i+ i + 5i 3ti + 5i itr -i-i) 
nii 



l(k + m r _4 + 5j 3 ,l + £/ 3 ,r-3 + ^4,2) 

m r _ 3 



(75) 

Here, the labels I2, h and I4 correspond to various factors in the coset SM ( 2 )i x ^( 2 )i x ^"( 2 )'--^ ^ namely the second factor 
su(2)i, the factor su(2) r _4 and su(2) r ^ 2 in the coset for l 2 , I3 and I4 respectively. These labels take the values 
l 2 = 1,2, Z3 = 1, 2, . . . , r — 3 and I4 = 1, 2, . . . , r — 1. We have written the functions y r in such a way that they are 
non-zero if (I2 + I3 + I4) mod 2 = 1. This fixes the suppressed label l\ corresponding to the first factor su(2)i in the 
coset to be l\ — (I2 +I3 + h) mod 2, due to the constraint h +l 2 + I3 = I4 mod 2. Finally, we introduced a 'generalized 
Kronecker-delta', 5 con d, which is 1 if the condition 'cond' is met, and zero otherwise. 

The matrix K = |A r _3 has rank r — 3 and is the Cartan matrix defined as above, the prime on the sum denotes 
the constraint that the summation variables are either even or odd, depending on the summation variable as well as 
the other labels of the functions y r . The parity of rrii is given by 



(k + I3 + 5i> r -l 4 5 r -2+h+i mod 2,0 + &i<l 3 5l 3 +i mod 2,l) Hiod 2 



(76) 



With this, we have completely specified the functions y r . To make the connection with the functions X ri 2k+i, there 
is one additional step, namely the introduction of the functions y r (k;l 2 ,l3,l4',q), which differ from y r (k;l2,l3,l4',q) 
only in the case that both I2J3 > 1. We define 



y r (k; 1, l 3 , l 4 ; q) = y r (k; 1, l 3 , 14, q) 

min(i 3 -l,2fc+l) 



(77) 



y r (k;2,l3,l4;q) 



1=0 



( Kfe+1 )/ 2 -([=-l )/4- r(fc _ m , 2 _ {l mQd 2) h _ ^ k . g) 



where \x~\ is the ceiling function, i.e. gives the smallest integer bigger or equal to x. The structure of the sum in the 
case that I2 = 2 is as follows. The argument V4 of the functions y r {k'; l 2 , 13, 14', q) is the same in all the terms. The 
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argument 1' 3 decreases in steps of one starting from 1' 3 = I3. The argument 1' 2 alternates between 1' 2 = 2 and 1' 2 = 1. 
Finally, the argument k' decreases by one every other term. We note that the value k' = — 1 can occur, in case that 
2k < l 3 - 2. 

With the functions y r in place, we can now write down the explicit expressions for the functions X r> 2k+i- We will 
do this for the ground states for u > first. There are four types of ground state patterns, namely (6, b — 1, 6, b + 1), 
(&, 6+1,6 + 2,6+1) and two others which are related by changing the heights from I to r — I, i.e. they are of the form 
(6, b + 1, b, b — 1) and (6, b — 1, b — 2, b — 1). Because the functions A" 2 fc+i remain unchanged when all the heights are 
reflected, we will concentrate on the first two ground state patterns for u > 0. In particular, we have the following 
result 

X rt 2k + i(a; b, b - 1, b, b + 1) = y r (k; 1,6-1, a; q) (78) 
X P)3 k+i(a; 6,6 +1,6 + 2,6 + 1) = y r (k;2,b,a;q) , (79) 

where we note that the functions ?/ r in the second line are a sum of terms y r in the case that 6 > 2. 

The ground state patterns for u < are of the form (6, 6 + 1, 6, 6 + 1) or (6, 6 — 1, 6, 6 — 1). By using the reflection 



symmetry (26 1 of the X r 2 k+i, we can restrict ourselves to the first set of patterns. As we explained in section 



HE 



the functions X r , 2 k+i(a; 6, 6 + 1, 6, 6 + 1; q) are related to the ones for the ground state patterns for u > 0, equations 



(32 ),(33 ). Namely, we have to distinguish two cases, the functions X ry2 k+i{a; 6+ 1, 6, 6+ 1, 6), where 6 = 1, 2, . . . , r — 3. 
The remaining case b = r — 2, i.e. X rt2 k+\{a] r — 1, r — 2, r — 1, r — 2) will be treated separately 

^r,2fc+i(a;6+ 1,6,6 + 1,6; g) = q^y r {k; l,6,a;g) (80) 
X r>2 fc+i(a; r - 1, r - 2, r - 1, r - 2; g) = q k+1 y r (k; 2,1, r - a; g) . (81) 

Finally, we come to the functions X r 2 k+i which are not related to ground state patterns, namely X r ^ 2 ], + i{a; b, b + 
1,6 + 2,6 + 3) and X r ^ 2 k+i{a; 6 + 3,6 + 2,6+1,6), and concentrate on the first one, where 6 = 1, 2, . . . , r — 4. As above, 



we use a relation to a ground state pattern for u > 0, namely (31 1, to obtain 



X rak +i(a; 6, 6 + 1, 6 + 2, 6 + 3) = q—y r {k; 2, 6, a; g) . (82) 

We numerically checked our conjecture for r = 5, . . . , 13, for all different ground state patterns up to sizes of at 
least k = 10 (in the case r = 13). We believe that it is possible to formalize the structure of the proof presented here 
for the cases r = 5,6, along the lines [3S], to prove the general case. 

VI. HEIGHT PROBABILITIES IN THE THERMODYNAMIC LIMIT 

In this section, we take the thermodynamic limit of the obtained form for the local height probabilities, in the case 



u > as well as for u < 0. In the former case, we can take the limit from the expressions ( 75 ) directly, while in the 
latter case, we first have to 'invert' the expressions, in order to obtain the result. For u < 0, the result allows us 
to interpret the local height probabilities in terms of CFT characters immediately. In the case u > 0, this is more 
complicated, and is the subject of the next section. 

A. Regime III (u > 0) 

We start by considering the case u > 0, and take the thermodynamic limit k —> 00, of the expressions for the height 



probabilities A ri2 fc+i(a; 6, c, d, e; g), which were expressed in terms of the functions y r (k; I2, h, Ia] q), given in eq. (78|. 
In those cases where both l 2 , l 3 > 1, these expressions are a sum over various y r {k'; 1' 2 , 13, 14,; q), but all terms, except for 
the first term y r (k] l 2 , l 3 , Z4; g), have an overall factor g zfe / 2 ; where I is some positive integer. So in the thermodynamic 
limit, we only have to consider the first term y r (k; l 2 ,l 3 ,l4;q), which is given in eq. (|75| 



Taking the limit lim^^oo y r (k; l 2 , l 3 , 1 4) q) is rather straightforward, because the k dependence only resides in the 
g-binomials, and we can make use of the result lim„_ 



j^— . The only complication lies in the parity of the 



n 
m 

summation variables, which depends on the parity of fc, forcing us to take the limit k — > 00 over either the even or 
odd integers. In the end, we obtain the result 



^lm-K-ra-lj 1<l4<r _ira r _i_i 4 



lim y r (k;l2,h,U;q)= J2 (aTla) II 

v WM1WK-3 i=2 



|(mi_i + m i+ i + 5i 3 a + <5z 4 , r -i-i) 
mi 



(83) 
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where the limit is taken either over even or odd k, and the prime on the sum denotes the constraints on the parities 
of the rrii, namely 

TUi = (k + l 3 + <5 J ;> I ._; 4 (5 r _2+/4+i mod 2,0 + <>i<l 3 Sl 3+ i mod 2,l) mod 2 . (84) 

The connection of these expressions with conformal field theory characters will be given in the next section. 

B. Regime II (u < 0) 



In order to obtain the thermodynamic limit of the local height probabilities for u < 0, we have to 'invert' the 
characters. The reason is that in the regime u > 0, the function </>(l), see eq. ( |18[ ), has to be maximized to find 
the ground states. Thus, to obtain the behavior at the critical point, we need to do the following procedure. For 
any finite system size k, the maximum value of 0(1) which can be obtained is ^(k + l)(k + 2). This means that the 
functions X 2 k+i we are interested in take the form q( k+1 ^ k+2 ^ 2 X 2 k+i(a; b, c, d, e; g" 1 ), see also [5 for more details. 



The particular boundary conditions relevant for the regime u < were described in section II D Equations ( 32 ) 



and (33 1 express the functions X 2 k+i for u < in terms of functions X 2 k+i relevant for u > 0. For the latter 



functions, we obtained the explicit fermionic expressions y r {k]l 2 ,l 3l l^q), given in equations (75) and ( |78| , but we 
note that for u < 0, the relevant functions have either 1% = 1 or 1$ = 1, so that we only need to consider the functions 



yr{k;l 2 ,kj4]q) in eq. (75 1 



For ease of notation, we write the functions y r {k]l 2 ^l 3 ,l^q) in the following way 

r-3 

,±m K m-iB m J~J 

i=l 



m;>0 
i=X,.,,r—3 



((1 



K)-m). + i(k + u) i 

771, 



(85) 



where the elements of the vectors k, u and B are given by ki = k(Si.i + <5i, r -3), = <5j,r-3<5z 2 ,i + <5i,/ 3 + <5i :1 — i— z 4 and 
Bi = Si r-i-u- The constraints on the sum are the same as those in (84). 



We are interested in calculating q^ k+1 ^ k+2 ^ 2 y r (k; l 2 , 13J4,; q 1 ), which is done by making use of the identity 



n + m 




n + m 






q —nm 






m 




m 











(86) 



and by changing the summation variables from the rrii to = | (k + u)i — (K • m)^. Most of the steps of this rewriting 
are straightforward. First, we note that we can write the rrii in terms of the rii as follows m 8 = (K _1 • (^k+ |u— n)) i; 
which allows one to write the function q( k+1 ^ k+2 ^ 2 y r (k;l 2 ,l3,l4 1 q~ x ) in terms of a sum over the rii. One has 
to take care in taking into account the constraints in the sum over the m, in the original expression. To this 
end, we consider the sum X^=i * n »> by making use of the explicit form of the matrix K, whose elements read 
Kjj = Sij — J5u_j'M, where i,j = 1, . . . ,r — 3. In particular, one finds Yn=i ' m )« = z ~§^ m r-3) giving rise to 

Ei=i in i = |(( r ~ 2 )( fc ~ nV-3) + EI=1 iu i)- 

The maximum value of irij is obtained for the minimal value r?i r _3, which is the parity of m r _3, which we denote 
as p r -3, i.e. Pr-3 = {k + I3 + $z 4 >3<5(z 4 -i) mod 2,0) mod 2. Thus, we find the following constraints 



r-3 



1 



i=l 
r-3 



r-3 

2)(k~p r - 3 ) + y^juj) 

i=l 

r-3 

2) (A: — p r -s) + iui) mod r — 2 



(87) 



where the second constraint follows from the fact that the parity of rn r -3 is fixed, which means that irii can only 
go down in steps of r — 2. 

Expressing the remainder of the expression in terms of the rii is straightforward, and leads to the following expression 

1/1 , 1W1 ,1 1 (*>+l)!2 X I r-l-21 3 ^ l 3 (r— 2-i 3 ) ( U - 1 ) ( r - 1 - 1 4 ) . 



r-3 



y^ v (k;l 2 ,h,k;Q) 



ni>0 



^n-K-^n-iB K- 1 !! J | 



((1-K- 



•n| 



i(K _1 • u)i 



(89) 
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where the prime denotes the constraints (87 1, and we note that the term — • K 1 ■ n can be written as 
— 5#2<; 4 <r-2(K _1 ■ n) r _i_; 4 . Finally the terms k in the binomials follow from |(K _1 ■ k), = k, for all i. The 
fact that all binomials contain the size dependence k is of great importance. It means that in the 'inverted' character, 
all particles are real particles, no pseudo particles are present anymore. Thus, in taking the thermodynamic limit 
k — > oo, all g-binomials are transformed in — . In this limit, the functions X 2 k+i(a; b, c, d, e; q) relevant for the 
critical behavior in the regime u < correspond to the characters of the primary fields of the Z r _2-parafcrmion CFT. 



To see this, we note that the first constraint in eq. (87) disappears in the limit k — > oo. The second constraint 



requires a bit more work. We first note that it is in fact independent of the size. Namely, k has the same parity as 
p r -3 if (I3 + <5z 4 >3<5z 4 -i m od 2,0) mod 2 = 0, which means we can write 

E, 1 N oN J ^ if (h + 6i 4 >3<5/ 4 mod 2,1) mod 2 = 1 

i(n t - - Ul ) mod (r - 2) = ■> 2 , . - 
,2 otherwise 

To simplify, we use the result that J2i^ u i = '3 + ( r ~ 1 — h)S2<u<r-2 + (r — 3)<5; 2 .i. The range on 52<i i < r -2 can 
trivially be extended to <$2<z 4 <r-i- The 'exceptional case' Z4 = 1 in the expression for im precisely allows us to 
extent the range ^2<z 4 <r-i to ^i<i 4 <r-ii i-e. the full range of Z4, by making use of the relation (l 2 + £3 + ^4) mod 2 = 1. 
Collecting all the terms, we finally obtain the simple expression for the constraint 



E 



l2+h~k-l 



Hi, 



mod (r - 2) . (90) 



We can now make the connection between the functions X r ^k+i in the limit k — > 00 and the characters of the primary 
fields of the Z r _ 2 parafermion CFT. For details on this theory, we refer to [H]. 
For b = 1, . . . , r — 3, we obtain the following result 

lim g3(*+i>(*+2)x 2fc+1 (a;&+ 1,6,6+1, 6; q' 1 ) = lim q^ k+1 ^ k + 2 ') y r (k; 1, b, a; q~ l ) = 



k- 



(91) 



(92) 



lim q 4 (— 2 > q 4 <*- 2 > yj. nv (fc; 1, 6, a; q) cx ch q ^_l 

k— >oo 

For 6 = r — 2, we obtain the following expression 

lim q^ k+1){k+2) X 2k+1 {a; r - 1, r - 2, r - 1, r - 2; q' 1 ) = lim q^ k+1){k+2) y r {k; 2, 1, r - a; q' 1 ) = 

k-^-oo k—too 

lim q a ^--^ y^ v (k] 2, 1, r - a; q) cx c\^ a r Z\ ■ 

k— >oo 

In the equations above, we denoted the fields of the Z,._2 parafermion CFT by & m , and the characters by ch q & m . 
These characters are given by the following expressions (see, for instance [TBH2TJ] ) 

<***J*= E s > (93) 

i—l,...r~3 

where the prime denotes the constraint X)I=i = m0Cl ( r — 2). We find that in the limit k — > 00, we obtain 
all the Z r _2 characters up to an overall factor of q. Namely, the labels l,m of the fields & m satisfy the constraint 
(I + to) mod 2 = 0. The label to is taken to be modulo 2(r — 2), because the fields & m = & m+2 ( r -2) are identified. In 
addition, the fields & m = &l^ 2 ~l 2 are also identified [44]. This implies that we can restrict the labels / and m to the 
range I = 0, 1, . . . , r — 2 and m = 0, 1, . . . , r — 3 for a total of |(r — 2)(r — 1) fields. Thus identifications above indeed 
cover all the fields in the Z r ^ 2 parafermion theory, establishing the connection between the latter and the critical 
behavior of the composite height model for u < (see also [5], noting that ch q & m = ch g $L m ). 

Before we continue with the details of the connection between the functions X 2 k+i for u > and CFT characters, 
we would like to point out the following. In inverting the characters above, we found that the pseudo-particles present 
in the fermionic description of the functions X r ^ 2 k+i become real. The reason behind this was that upon inversion, all 
the q-binomials acquired dependence on the size fc, and therefore became factors of l/(q) n in the limit k — > 00. This 
is in fact the generic behavior. In the original ABF model, the critical behavior for u > is given in terms of minimal 
models, which have a fermionic description, with one real particle, while the other particles are pseudo-particles. 
The matrix K was the same as the one used here. Upon inversion, all the pseudo-particles become real. In fact, 
many (anyonic) chain models exhibit regions whose critical behavior is described in terms of minimal models. Upon 
inverting the sign of the hamiltonian, one often finds another critical region, whose criticality is described in terms of 
the Z r _ 2 -parafermion theory. One can speculate that this behavior is governed by integrable points, which exhibit 
the same behavior as observed for the ABF model, and the composite height model considered here. 
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VII. CONNECTION WITH A COSET CFT: REGIME III (u > 0) 



In this section, we explain the connection between the expressions for the local height probabilities in regime III 
(u > 0) and characters of a particular coset conformal field theory. This connection for regime II (u < 0) was explained 
in the previous section. 

To make the connection between the explicit UCPFs which we obtained for the functions X Ty 2k+i, the local height 
probabilities, and the CFT describing the critical behaviour of the model, we consider bosonic forms of the characters 
associated with the coset conformal field theory "^'"^j'^^'"' ■ 

Before we start, we first collect a few properties of this coset theory. The fields $ of this theory carry, a priori, 
four labels, associated with the four su(2) algebras. These four labels, taking the values l\, 1% = 1, 2, Is = 1, . . . , r — 3 
and I4, = 1, . . . , r — 1 satisfy the constraint l\ + I2 + ^3 + I4 = mod 2. Throughout, we make the choice l± — 1, and 
frequently omit this label, and write the fields as ^ > ^'' 3 - 

The scaling dimensions ft.i 2 ,z 3 ,z 4 of the fields were obtained in [3], making use of the Coulomb gas results of [55], in 
particular 

»mm = ,,,4;: y 4+ h+ L>.i' (94) 



8r(r-2) ' 8 



The scaling dimensions satisfy /i(3 — Z2, r — 2 — Z3, r — li) — h{l%, h, h), reflecting the fact that the fields $^-24 3 l2 ' r 2 ' 3 

and (J)? 1 '' 2 '' 3 are identified. 

4 

One way to view this coset theory is via a product of two unitary minimal models. The unitary minimal models 
have a coset description, su ^ 2 l^2"^ r ~ 3 being the minimal model Ai(r — l,r), where r — 4 corresponds to the 
Ising CFT. The coset SM ( 2 )i x ^^j 1 x s"( 2 )r-4 can ^ e v j ewec j as ^j^g product of two minimal models, in particular 
M(r - 2,r - 1) k M(r - 1, r) = s " (2 ^^ (2 3 } - 4 k su(2) l fi 2 ™™ r - 3 ■ The coset g^gkjlggli^g^ is not the direct 
product of the two consecutive minimal models, but corresponds to a non-diagonal modular invariant of the product 
theory. We note that in [45], the case r = 5 was considered in the context of non-abelian quantum Hall states. 

To construct this modular invariant, we will use the intuition that one can construct 'new' CFTs from 'old' ones 
by condensing a boson (i.e., a particle corresponding to a field with integer scaling dimension) present in the theory, 
as advocated in |46j . Under this condensation, the boson is identified with the vacuum, or identity field; in addition, 
those fields which are related to each other by fusion with the boson (i.e., the condensate), are to be identified. Two 
other steps are necessary to construct the new CFT. First, two fields which are each others dual (i.e., their fusion 
contains the identity) have to be split, if both the identity and the boson which is condensed are present in the fusion 
product of the two fields considered. This is because the boson is identified with the identity, and two fields can not 
be fused to the identity in more than one way. Second, fields with have non-trivial monodromy with the boson, are 
confined. We used these principles to guide us in constructing a non-diagonal modular invariant of the product theory 
of two consecutive minimal models. More details about the condensation procedure can be found in [4*5] . 

To construct the relevant modular invariant, we first briefly recall some facts about minimal models lj. Minimal 
models are labeled by two co-prime integers r and r', where we choose the ordering r < r' . The minimal model 
M(r,r') is unitary if r' = r + 1. The fields (f> present in this theory can be labeled by two integers (£3, £2); where 

1 < h < r ' — 1 an d 1 < I2 < r — 1. The scaling dimensions of the fields </>(z 3l ; 2 ) are given by hn 3 ,i 2 ) = — 4rr' {r r ' . 
One finds that hn St i 2 \ = hr r '-i 3! r-i 2 )i and the labels (hjh) and (/ — l^,r — 1%) indeed correspond to the same field, 
1 (h,h) — 0(r'-J 3 ,T— 2a)- Thus, the number of fields in the model M(r,r') is |(r — l)(r' — I). 



The chiral characters of the minimal models A4(r,r') can be written as 



* £li (<Z) = £ q n(nrr'+l 3 r-l 2 r>) _ ? (nr+2 2 )(nr' +l 3 ) ( (95) 



where c is the central charge of the corresponding theory, and h the scaling dimension of the field. In addition, 

(?)oo = n~=i(i -<?")■ 

The fields in the product theory M. (r — 2, r — 1) X M. (r — 1, r) are, in anticipation of the results given below, labelled 
by (?2) £3; h, '2)1 and the scaling dimensions are given by the sum of the scaling dimensions of the fields in the two 
factors. In particular, the field (3, 1; 1,3) has scaling dimension /i(3,i ; i,3) = 2, irrespective of the value of r (as long 
as r > 5, which is necessary for the coset to be defined). By applying the condensation strategy outlined above, one 
can find a set of fields which are not confined, and are inequivalent of one another. This set of fields has the labels 
(1, Is; 2*4, 1) in the case that Is mod 2 = 1 and (2, 2*3; I4, 2) when Is mod 2 = 0. We note that Is and I4 take the values 
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Is = 1,2, ... ,r — 3 and I4 = 1, 2, . . . , r — 1, which means that the product theory contains (r — 3)(r — 1) different 
fields. The fields in the original product theory which are identified with the field with the label (1, Z3; I4, 1) are of the 
form (I, Z3; I4, 1), with I mod 2 = 1. The fields identified with (2, l 3 ; I4, 2) also take the form (I, l 3 ; I4, 1), but now with 
I mod 2 = 0. 

We will denote the characters of the field in the product theory by X(; 2 ,z 3 ;Z 4 , ;;,)(<?) = X(i 2 ,h) (l)X(u,i' 2 ) (<?)• With these 
fields, one can construct different modular invariants. As usual, there is the diagonal invariant, corresponding to the 
direct product of the two minimal models 

(r-2,r-3:7'-l,r--2) 

Zdiag = \X(i 2 ,i 3 M' 2 )(<l)\ 2 1 ( 96 ) 

(!a,J 3 ;J4,«i)=(i,i;i.i) 

where the sum is over all the fields in the product theory. 

Apart from this diagonal invariant, one can construct a different modular invariant, which corresponds to the coset 
sn(2)ixsn(2)ixs-u(2) r _ 4 ^ j en0 ^- e ^is theory by Ai(r — 2, r — 1, r), and the partition function, which was obtained 
from the condensation picture outlined above, takes the form 

r— 3 1 — 1 r— 3 r— 1 

Zm{t-2,t-X,t) = E IX(l.W4, 1)0?) + X(3,i 3 ;Z 4 ,3) OH" '| 2 + E E IX(2,i 3 ;i 4 ,2) (?) + X(4,i 3 ;i 4 ,4) (?) + ' ' ' ? • (97) 

i 3 = l Z 4 = l Z 3 =2 Z 4 =l 

I3 odd Z 3 oven 

This form of the modular invariant gives expressions for the characters of the coset M.(r — 2,r — 1, r), in terms of the 
bosonic characters of the minimal models. The characters corresponding to the fields in the coset model can now be 
written in terms of the bosonic characters. 

The characters of the coset fields are denoted as ch g $^'' 3 and are given by 

r-2 

c\$£' h = E X(W,;U,o(«), ( 98 ) 
1=1 

l+l 3 =0 mod 2 

where it is assumed that (I2 + I3 + I a) mod 2 = 1. We checked numerically that these characters indeed correspond to 
the branching functions of the coset M.{r — 2, r — 1, r). 

We can now relate the thermodynamic limit of the functions X r ,2k+i to the characters of the coset fields. As 
already indicated in [5], one has to consider the cases k odd and k even separately. The reason is that the parity 
of the summation variables in the functions y r {k;l2,lz,l4',q) depends on the parity of k. We obtain the following 
identification. For k — 2p + 1 odd, the limits are, for the ground state patterns of type and Gf , 

lim X rAp+3 (a; b + 1, 6, b + 1, 6 + 2; q) = lim y r {2p + 1; 1, b, a; q) = ch„$^ h (99) 

p— >oo p— too 

lim X rAp+3 {a; 6, b + 1, b + 2, b + 1; q) = lim y r (2p + 1; 2, b, a; q) = ch„$^6 ( 10 o) 

p— >oo p—toc 

= lim y r (2p + 1; 2, r - 2 - b, r - a; q) = ch^^ 2 ~ b . 

p— f 00 

For k = 2p even, the connection between the functions X r _2k+i for the u > ground state patterns is slightly 
different, namely 

lim X rAp+1 {a;b+l,b,b+l,b + 2;q) = lim y r (2p; 1, r - 2 - b, r - a; q) = di^l^ 2 ^ (101) 
lim X r , 4p+ i(a; 6,6 + 1, b + 2,6 + 1; g) = Km y r (2p; 2, 6, a; g) = ch q $ 2 > b (102) 

= lim 2/ r (2p; 2,r - 2 - b,r - a;q) = ch q <S> 2 r ^ 2 ~ b . 

p—too 



In taking the limit k — > 00, the structure of the resulting expressions for y r is given in Eq. (83 1. The first and 
last q-binomials are transformed into l/(q) mi and l/(q) mr _ a . With this, the equivalence between the two different 
identifications for the ground state patterns Gf follows from a simple re-parametrization of the sum over the mj. 
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VIII. DISCUSSION 



While the anyonic chains were introduced as a simple setting to study interacting anyons appearing in various 
topological phases, it has turned out that they have rich phase diagrams interesting in their own right, much like 
for ordinary spin chains. The two integrable critical points identified in Ref. [5] led to the composite height model 
we have been studying here. In addition to novel phase diagrams and critical behavior, to mention one thing, the 
one-dimensional anyonic chains are of direct relevance to the boundary behavior of a nucleated phase, arising from the 
interactions, and the 'parent' topological phase hosting the bare anyons [IS]. By now, transitions between different 
topological phases due to the interactions between the anyonic excitations in topological liquids has been studied in 
quite some detail, see for instance 45, 48, 49 . 

In this paper, we have further investigated the properties of the composite height model of Ref. [5J as follows. We 
obtained and studied the fermionic forms of the LHPs and identified their off-critical CFT structure at two different 
regimes. The same CFTs are related to the critical points of the height model and that of the original anyon chain. 
The proof was based on the recurrence properties of the LHPs and the UPCFs, very much like in the original case of 
the ABF model and the minimal models. Proofs of this type, based on polynomial recurrences, are straightforward 
and tractable but unfortunately much of the physics, especially the off-critical CFT structure, in the height model is 
obscured as a trade off. 

In particular, we gave a closed UCPF form for the LHPs. Using these, we were able to analytically prove the 
correspondence for r — 5,6 and gave a general conjecture based on the structure of the UCPFs for r — 5,6, correct 
asymptotic central charges, as well as numerical checks. Although the form of the recursion we used here quickly 
becomes cumbersome as r grows larger, we suspect that our proof can be further improved to a proof for general r, 
along the lines of [25], using the UCPFs conjectured here. We note that the fermionic representations of UCPFs as 
unitized characters are usually not unique [13], as there might be several integrable perturbations of the CFT and 
different ways to introduce the 'finitization' in the size k. We have not attempted to analyze to what integrable 
perturbations our UPCFs correspond to. Moreover, the recursions we were forced to use are more general than those 
in the lattice model, and the physical quantities related to the height probabilities were only obtained as special cases 
and, in the regime III, as linear combinations from the functions y r =5,e- I n addition, we obtained various relations 
between the general functions, which were necessary to show the equivalence with the local height probabilities. We 
note that in the case of the original ABF models, one can show more directly that the recursions for the local height 
probabilities and the UCPFs are identical. 

To the best of our knowledge, the fermionic forms for the characters of the coset su ( 2 )i xsu (?h xsM(2) r _ 4 a 
and characterized by the fact that they have two real fermions. In addition, the forms based on the characters of 
minimal models also appears to be new. From these characters, we obtained characters of the Z r _2 parafermions by 
sending q — > q^ 1 and taking the limit of large systems size in the obtained UPCFs. In the context of the anyonic 
chains, this procedure corresponds to changing the overall sign of the hamiltonian. These infinite system size Z r _2 
parafermion characters are of course not new, but the unitized versions we obtained here do differ from the unitized 
Z r _2 parafermion characters which correspond to regime II in the original ABF model. The latter can be obtained 
from the finitized characters corresponding to regime III of the ABF model by the inversion procedure, see for instance 
[25] . It is interesting that the finitized characters of the diagonal coset models generically lead to different finitized 
characters of the Z r _i parafermions! 

We briefly remark on the modular properties of the fermionic character formulae for the coset theories presented 
here. As our analysis in section [V] shows, a crucial role is played by the pseudo-particles |37j . In the limit k — > oo for 
u > 0, our fermionic formulae are generalized r — 3-dimensional g-hypergeometric series containing r — 5 g-binomial 
factors with finite arguments from the pseudo-particle sector. In the case when there are no pseudo-particles, Nahm 
|50j has provided a conjecture regarding the modular properties of multi-dimensional g-hypergeometric series arising 
from UCPFs determined by a bilinear form K and a 'shift' B. In the case r — 5, there are no pseudo-particles present, 
and our formulae give back the well-known results and (asymptotic) central charges c = |, | based on K = ^A 2 , ^A^ 1 
at rank 2 [5] [50]. In the general case for u > 0, our fermionic character formulae (83) have g-binomial contributions 
from the pseudo-particles, and the modular properties of the resulting g-series are even more complicated than in 
Nahm's conjecture and left for future study. 

Lastly, we would like to mention some future directions for the study of the composite height model that were left 
out from this paper. Much like in the original ABF models, one would like to obtain bosonic forms for the characters 
of the coset theory, to pave way for corresponding Rogers-Ramanujan type identities and the modular properties of 
the coset theory. The UCPFs and their recursions that we have explicitly used in our proof for r — 5,6 are more 
general than the LHPs and characters, that are obtained only at special values of the arguments and, for the coset 
theory as linear combinations. This behavior is new compared to the ABF models and it would be interesting to 
find the natural representation theoretic setting, if any, of these q-polynomials or, conversely, to find a more direct 
functional form for the recursions in the LHPs. To study the physics and combinatorics of the height model more 



29 



directly, a study of LHPs in terms of lattice paths of the composite height model [5TJ [32], where (181 acts as the 
Virasoro generator Lq, would be interesting. This would at the very least combinatorially relate generating functions 
in the path space of the height model to our fermionic UCPFs. Also, we have addressed only a half of the phase 
diagram of the composite height model, < p < 1, along the lines of the original paper [5], that corresponds to the 
quantum mechanical anyon chain. The negative p regime, —1 < p < 0, can also be studied via the CTM method but 
the quantum mechanical interpretation, if any, is unknown. The details of the phases of the composite height model 
for negative p will be given in a subsequent publication [53] . 
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Appendix A: Recursions for r — 6 



Here we present the proof of the equivalence of the remaining recursions for r — 6, one for each a = Z4 for the different 
patterns determined by the height b. The other values for a = I4 follow in almost the same way, since the patterns in 
the recursions and the indices l 2 ,l 3 are the same irrespective of a = I4, and the identities for y( ABC \k\,k 2 ,k 3 ]l 2 ,l 3 ,l4) 
are identical for fixed parity of a = Z4; thus the recursions differ only slightly in the specific parities appearing. 



1. Recursions in Gf 

The recursion for X2fc+i(2; 1232) is 

X 2k+1 {2: 1232) = X 2fe _ 1 (2; 1212) + X 2fe _ 1 (2; 3212) 



which is 



Now 



X 2k+1 (2; 1232) = q k y( BBB \k - 1, k - 1,0; 212) + y {BBB \k - 1, k - 1,0; 234) 
+qiy( AAB )(k - 2, k - 2, 0; 124) + q^y^ ABA \k - 2, k - 2, 0; 214). 



i/ ABC) (h,k 2 ,k 3l 234) = y( CBA \k 2 ,k l7 k 3 ;212) 
y (ABC) (ki,k 2l k 3 ;124) = y^ CBA \k 2 , k u k 3 + 1; 212) 
y {ABC) (ki,k 2l k 3 ;2U) = y^ CBA \k 2 - 1, fc x + 1, fc 3 ; 212). 



These give 



X 2k+l {2; 1232) = q k y^ BBB \k - 1, k - 1, 0; 212) + y {BBB \k - 1, k - 1, 0; 212) 
+qiy( BAA )(k - 2, k - 2, 1; 212) + q^y^ ABA \k - 3, k - 1, 0; 212). 

We use the relation in ( [63] ), 

2/ (ABC) (fci, k 2 , k 3 ; 212) = y ^+ 1B + 1A+1 \k 2 -l,h + 1, k 3 ; 212), for k 3 > 0, 

to get 

X 2k+1 (2; 1232) = q k y^ BBB \k - 1, k - 1,0; 212) + y^ BBB \k - 1, k - 1,0; 212) 
+qh {BBA) (k - 3, k - 1, 1; 212) + q^y^^ik - 3, k - 1, 0; 212). 

Now apply the recursion 

y(BBA) _ Sk _ ll . 212 ) = y(BBA) q, _ 3; fc _ l5 _ 1; 212) + y^ BAA \k - 2, k, -1; 212) 



(Al) 
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to get 

A 2fe+1 (2; 1232) - q* {y^ BBA \k - 3, fc - 1,-1; 212) + q^y^ ABA \k - 3, k - 1,0; 212)) 
+q k y (BBB) (jfe - l, jfe - 1, 0; 212) + y (BBB \k - l,fe - 1,0; 212) + qiy {BAA ^ (k - 2, jfe, -1; 212). 
This is equal to 

A 2fc+1 (2; 1232) = qiy {BBA \k - 1, * - 1, -1; 212) + <zV SBS )(fc -1,^-1,0; 212) 
+y( SBS )(fc - 1, k - 1, 0; 212) + qiy^ BAA \k - 2, jfe, -1; 212). 

The RHS is equal to 

g l (y<- BBA ) (jfe - 1, fc + 1, -1; 212) + (jfe - 2, jfe, -1; 212)) + y^ 8 ^ [k - 1, k - 1, 0; 212) 

= q V 5 - 4 ^ " 2, fe, 1; 212) + y^ AAA \k - 2, k, 0; 212), 



where we use ( 63 ) again on the second term. So we are left with 

X 2k+1 {2- 1232) = qiy^ BAA \k - 2, k, 1; 212) + y^ AAA \k - 2, k, 0; 212) 

which is the original recursion for X 2 fc+i(2; 1232) = y( AAA \k, k, 0; 212). The recursion for A 2 j, +1 (4; 1232) is similar 
and omitted. 

Finally, the recursion for X 2k +i{^] 3454) is 

X 2fc+1 (4; 3454) = q^X 2fe _ 1 (4; 1234) + q^ X 2k _ 1 (A; 3234) + X 2k _ 1 (A; 3434) + X 2fc _ 1 (4; 5434). (A2) 
This is the same as 

X 2fe+1 (4; 3454) = q^qiy^ BAB \k - 1, jfe - 1,0; 214) + q^y^ BBA \k - 1, k - 1, 0; 124) 
+qiy {BAA \k - 1, k - 1, 0; 122) + y {BBB \k - 1, jfe - 1,0; 212). 

Next, we use the relations on the last the two terms, which are not part of X 2 fe + i(4; 3454), 

y^ ABC Hk 1 ,k 2 ,k 3 ;122) = y^ 8 ^ (k 2 + 1, fc x - 1, fc 3 + 1; 234) 
y {ABC) {k u k 2 ,kz\ 212) = y (CSA) (fc 2 ,fc 1 ,A : 3;234). 

So we get 

X 2fe+1 (4; 3454) = g^g V^ B) (* _ i, fc - 1,0; 214) + q^y^ BA \k - 1, k - 1,0; 124) 

+g l y (^B) ( fc> jfe _ 2) 1; 234) + y( SBS ) (jfe - 1, jfe - 1, 0; 234) . 

Also, using the relation 

2/ (ABC) (fci, k 2 , k 3 ; 234) = y( c + 1B + 1A + 1 )(fe 2 + 1, * x - 1, fc 3 ; 234), for fc 3 > 0, 



in ( 63 ) on the last term, the recursion reduces to 



X 2fe+1 (4; 3454) = q^qiy^ BAB \k - 1, k - 1,0; 214) + q h ^y { ~ BBA \k - 1, k - 1,0; 124) 

+(? l 2/ (^B)( fc; jfe - 2, 1; 234) + (jfe, jfe - 2, 0; 234) 

which is simply the same as 

X 2fc+1 (4; 3454) = q^y {BAB \k - 1, jfe - 1, 0; 214) + q^y( BBA \k - 1, k - 1, 0; 124) + y (AAA) (k, jfe, 0; 234), 
as desired. Again, the recursion for X 2k+ i(2; 3454) is similar and left for the reader. 
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2. Recursions in G 2 

The recursions for X 2k+ i(2; 3234) is 

X 2k+1 (2; 3234) = X 2fe _ 1 (2; 1232) + q^ X 2k ^(2; 3232) + X 2fe _ 1 (2; 3432) + X 2k ^(2, 5432). (A3) 

Writing the RHS in terms of y's gives 

X 2k+1 (2; 3234) = q^y^ cc \k — 1, fc — 1, 0; 212) + q^ q^ y^ c AA \k - 1, fc - 1, 0; 122) 
+?/ c ' CM >(fc - 1, fc - 1, 0; 124) + qiy^ CAC \k - 1, fc - 1, 0; 214). 

Using 

y( ABC )(fc 1 ,fc 2 ,fc 3 ;212) = j/ ASC > (fc x + 1, fc 2 - 1, fc 3 - 1; 122) 
y^ ABC Hh,k 2 ,k 3 ;l2A) = y^ CBA \k 2 + 1, k x - 1, fc 3 ; 122) 
y( ABC )(fc 1 ,fc 2 ,fc 3 ;214) = y( CSA ) (fc 2 , fc l7 fc 3 - 1; 122), 

we get 

X 2k+1 {2- 3234) = q^y( ccc \k, k - 2,-1; 122) + y^ ACC \k, fc - 2, 0; 122) + qiy^ CAG \k - 1, fc - 1, -1; 122) 

+q^qi y(° AA ){k - 1, fc - 1,0; 122). 

Next we use the recursion on the last two terms 

q^y( CAA \k - 1, fc - 1, 0; 122) + y( CAC \k - 1, fc - 1,-1; 122) = y^ CAC \k - 1, fc + 1, -1; 122), 



so 



X 2fe+1 (2; 3234) = q^ (q^y^ ccc \k, k - 2,-1; 122) + y {CAC \k - 1, k + 1, -1; 122)) + y (ACC) (fc, fc - 2, 0; 122). 



Focusing on the hrst term, we have by ( 63 ) and the recursion 

qhyiCCC)^ k _ 2) _ 1; 122) = ? l/2 y (AAA) (fc _ 1? fc _ 1; _ 1; 122 ) = 
q l/2 y (AAA)^ + l k _ 1( _ 1; 122) _ q ^ y (CAA)( k _ l k _ 10 . 122) 

and 

y (CAC) (fc _ X) fe + ^ _ 1; 122) = y(ACA)^ k _ 2j _ 1; 122) + ^yfcci)^ fc _ 2 , 0; 122). 



Using ( 63 ) again, two terms cancel and give 

q i/2 y (AAA) (fc + 1, fc - 1, -1; 122) + y( ACA ^ (fc, fc - 2, -1; 122) = y^ ACA ^ (fc, fc - 2, 1; 122). 

Finally we get 

X 2fe+1 (2; 3234) = qiy ACA {k, fc - 2, 1; 122) + y (ACC) (fc, fc - 2, 0; 122) . 

The recursion for X 2k+ i(A; 3234) is similar and omitted. 
Finally, the recursion for X 2 fc + i(3; 4345) is 

X 2fe+1 (3; 4345) = q^ 1 X 2k ^(S; 2343) + q^ X 2k _ 1 (3; 4343)+ X 2fc _i(3; 4543) . (A4) 
In terms of y's the RHS is 

A 2fc+1 (3; 4345) = q^ (y^ BAB ^ (fc - 1, fc - 1, 0; 223) + qiy^ AAB ) (fc - 2, fc - 2, 0; 113)) 
+q h ¥qi y i BAA ) (fc - 1, fc - 1, 0; 133) + y {BBA ^ (fc - 1, fc - 1, 0; 113). 



32 



Now 

y^ ABC \k 1 ,k 2 ,k 3 ;223) = y^ ABC \k u k 2 - 2, k 3 + 1; 133) 
y iABC \h,k 2 ,k 3 ; 113) - y (ASC) (fc 1 + l ! fc 2 -l,fc 3 ;133) 

and we get 

A 2fe+1 (3; 4345) = q^y {BAB) {k - 1, k - 3, 1; 133) + V AAS) (fc - 1, fc - 3, 0; 133) 

+(Z i ^ (Z l 2; (^A) (fc _ 1? j. _ 1(0 . 133) + y ( BS ^)(fc, fc - 2, 0; 133) . 

Now, using the identity y (ABC) (h, k 2 , fc 3 ; 133) = y {CBA) {k 2 + 2, k x - 2, fc 3 ; 133) in ([59]) on the last term, gives 

X 2fe+1 (3; 4345) = y {ABB) {k, fc - 2, 0; 133) 
(y (BAB) (k - 1, k - 3, 1; 133) + g V AAS) (fc - 1, fc - 3, 0; 133) + q^y {BAA \k - 1, fc - 1, 0; 1133)) . (A5) 

The term in the brackets is, using the recursion in the first term and ( |59[ ) in the second, 

y(BAB) (fc _ 1} k _ 3 1; 133) + q i y (AAB)^ _ lk _^Q. 133 ) + ( fc _ ^ fc _ ^ Q. H 33 ) = 



y 



{BAB) 



(fc - 1, fc - 3,-1; 133) + y {BBB) (k, k - 2, -1; 133) + qiy {BAA \k - 1, fc - 3, 0; 133) 
+<z !j/ s ^)(fc- l,fc- 1,0; 133), 



or further, using eqs. (|59j) and (62 1 on the last term, 

y {BAB ){k _ lk _ ^ _ 1; 133) + y(BBB)^ k _ ^ _ 1; 133) + q ^ y (BAA)^ _ 1} jfe _ X> Q; 133) 

= y( SAB )(fc - 1, fc - 1, -1; 133) + y {BBB) {k, k - 2, -1; 133) + qiy {AAB \k - 1, fc - 1, 0; 133) 



(BAB) 



(fc + 1, fc - 1, -1; 133) + y (BBB) (k, fc - 2, -1; 133) = y (BBB) (k, fc - 2, 1; 133) 



= j/( ASA )(fc,fc-2,l;133) , 

which follows again by (62). Returning to (A5), we get back y( ABB \k,k,0;133) = A 2fe+ i(3; 4345) as desired. The 
recursions for X 2 fc + i(l; 4345) and X 2 j.+i(5; 4345) are similar and omitted. 
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